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FOREWORD 


Rocketdyne,  a  Division  of  North  American  Aviation, 
Inc.,  has  prepared  this  final  report  which  documents 
in  detail  all  of  the  work  performed  in  fulfillment 
of  the  program  “Boiling  Heat  Transfer  for  Cryogenics” 
sponsored  by  the  National  Aeronautics  and  Space 
Administration,  George  C.  Marshall  Space  Flight 
Center,  Huntsville,  Alabama,  under  Contract  No, 
NAS8-5337*  The  work  was  performed  by  the  Heat 
Transfer  and  Fluid  Dynamics  Unit  of  the  Research 
Department  during  a  9-nionth  period  extending  from 
l6  May  1963  to  15  February  1964. 
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ABSTRACT 


I '■  An  extensive  survey  of  available  information  from 
'•  ^  /f.V' 

the  literature  on  heat  transfer  to  boiling  hydrogen, 

nitrogen,  and  oxygen  is  presented.  Included  are  a 

bibliography  of  available  experimental  data,  a  summary 

of  boiling  theory,  graphical  comparisons  of  the  experi- 

mental  data,  and  compilations  of  pertinent  physical 

properties  for  each  of  the  three  cryog^hic  fluids/ 


Based  upon  the  results  of  the  survey,  recommendations 
are  made  for  a  test  program  to  provide  critical  infor- 

f  '■ 

mation  where  data  are  lacking.  ' 
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SUNMAEY 


This  report  presents  an  extensive  survey  of  the  available  information 
from  the  literature  on  heat  transfer  to  boiling  hydrogen^  nitrogen,  and 
oxygen  for  both  pool  and  flow  conditions.  A  comprehensive  bibliography 
of  74  literature  articles,  reports,  unpublished  papers,  and  theses  con¬ 
cerning  experimental  data  is  presented.  Each  data  reference  is  coded 
with  respect  to  boiling  regime,  flow  condition,  geometry,  unusual  condi¬ 
tions,  and  method  of  data  presentation.  From  this  coding,  a  large  number 
of  convenient  cross-reference  bibliography  lists  were  prepared  for  spe¬ 
cific  conditions  of  boiling.  In  addition  to  the  coding,  each  data  ref¬ 
erence  is  rated  with  respect  to  the  experimental  technique  description. 

A  review  of  all  the  available  theoretical  approaches  to  the  prediction 
of  boiling  phenomena  is  presented.  In  general,  empirical  expressions 
are  not  considered.  Calculations  based  on  the  theories  are  given  graph¬ 
ically  for  each  of  the  three  cryogenic  fluids  for  several  types  of  boil¬ 
ing,  Computations  with  the  equations  were  made  possible  by  an  extensive 
compilation  of  all  the  physical  properties  pertinent  to  boiling. 

The  available  experimental  data  are  presented  in  graphical  form  and 
compared  to  theoretical  predictions  v/here  possible.  The  effects  of  the 
many  boiling  parameters  involved  are  discussed.  Based  on  the  results 
of  the  literature  survey,  recommendations  are  made  for  additional  ana¬ 
lytical  and  experimental  studies  that  are  required  to  provide  an  adequate 
knowledge  of  heat  transfer  to  boiling  cryogenics. 
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INTRODUCTION 


Interest  in  heat  transfer  to  boiling  cryogenic  fluids  such  as  hydrogen, 
nitrogen,  and  oxygen  has  increased  rapidly  during  the  past  15  years. 

This  interest  has  been  due  primarily  to  a  growing  number  of  space 
applications o  Compared  to  the  boiling  of  ordinary  fluids  such  as  water 
at  temperatures  near  or  above  ambient,  cryogenic  boiling  is  character¬ 
ized  by  smaller  temperature  driving  forces,  smaller  heat  fluxes,  and 
the  practicability  of  operation  in  the  stable  film  boiling  regime. 

This  is  illustrated  graphically  in  Fig.  1  where  the  atmospheric  pool 
boiling  curve  of  water  is  compared  to  hydrogen  and  oxygen  for  horizontal 
wires.  As  can  be  seen,  nucleate  boiling  of  hydrogen  occurs  for  A  T 
values  to  only  about  5  F  and  to  heat  fluxes  on  the  order  of  10,000 
Btu/hr-sq  ft.  For  oxygen,  the  nucleate  boiling  range  is  extended  to  a 
A  T  of  about  20  F  and  a  heat  flux  of  about  50,000  Btu/hr-sq  ft.  For 
water,  the  nucleate  boiling  A  T  is  greatly  extended  to  approximately 
50  F  and  the  heat  flux  to  approximately  400,000  Btu/hr-sq  ft.  The 
direct  transition  of  nucleate  boiling  to  stable  film  boiling  for 
hydrogen  and  ox^^gen  results  in  surface  temperatures  below  1000  F,  which 
is  attainable  with  many  materials.  Thus,  heat  fluxes  greater  than  the 
pool  nucleate— boiling  maximum  heat  flux  can  be  achieved  without  burnout. 
With  water ,  however ,  direct  transition  to  stable  film  boiling  can  lead 
to  surface  temperatures  in  excess  of  2000  F.  For  many  materials  this 
will  result  in  a  burnout,  and  operation  consequently  is  restricted  to 
the  nucleate  boiling  regime. 

These  characteristics  of  cryogenic  boiling  complicate  and  make  difficult 
the  acquisition  of  accurate  and  complete  experimental  data.  Therefore, 
it  seems  appropriate  to  survey  the  available  data  and  determine  the 
status  of  knowledge  of  cryogenic  boilingo 
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Heat  Flux,  q/A,  Btu/hr~ft 


FIGURE  1 


COMPARISON  OF  POOL  BOILING  CURVES 


This  report  presents  such  a  survey.  Included  is  a  complete  bibliography 
of  available  experimental  data.  The  references  are  coded  and  cross- 
referenced  to  facilitate  searches  for  specific  information.  Also 
included  are  compilations  of  pertinent  physical  properties,  a  summary 
of  available  boiling  theory,  and  graphical  presentations  of  the  available 
data.  Based  upon  the  results  of  the  survey,  recommendations  are  made 
for  a  test  program  to  provide  critical  information  for  regions  where 
data  are  lacking. 


LITERATURE  SURVEY 


An  exhaustive  search  of  the  literature  was  conducted  for  all  literature 
articles,  reports,  unpublished  papers  and  theses  concerning  experimental 
data  for  boiling  heat  transfer  from  solid  surfaces  to  hydrogen,  nitrogen, 
and  oxygen.  Most  of  the  available  literature  to  I96O  was  surveyed  by 
Richards,  Steward,  and  Jacobs  (Ref,20l).  The  present  survey  includes 
all  articles  to  January  1964,  The  sources  consulted  in  the  literature 
survey  were: 

Chemical  Abstracts  (1907-June  I962) 

Heat  Transfer 
Each  Compound 

Nuclear  Science  Abstracts  (1946— December  1963) 

Heat  Transfer 

Boiling 

Each  Compound 

Engineering  Index  (l944-February  1964) 

Heat  Transfer-Boiling  Liquids 

Applied  Science  and  Technology  Index  (1958-January  1964) 

Boiling 

Low  Temperature 
Heat  Transmission 
Each  Compoimd 

Industrial  Arts  Index  (1941-1957) 

Boiling 

Heat  Transmission 
Each  Compound 

Aeronautical  Engineering  Index  (1947-1956) 

Heat  Transfer 

Aerospace  Engineering  Index  (1957-1958) 

Heat  Transfer 
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International  Aerospace  Abstracts  ( 19 6 1- January  1964) 

Heat  Transfer 

Refrigeration  Abstracts  (1949-^957) 

Dissertation  Abstracts  (l995'“December  1963) 

Liquid  Propellant  Information  Agency  (195^-1962) 

Heat  Transfer 
Cryogenic  Liquids 
Each  Compound 

Chemical  Propulsion  Information  Agency  (1963) 

Armed  Services  Technical  Information  Agency  (ASTIA) 

Special  bibliography  vas  requested  on  Boiling  Heat  Transfer. 
The  bibliography  was  compiled  from  ASTL4  holdings  from  1952  to 
January  1964. 

NASA--NACA: 

NASA  Scientific  and  Teclmical  Aerospace  Reports  (STAR) 

January  1963  to  January  1964 

NASA  Technical  Announcement  Bulletin 
April~Dec ember  1962 

NASA  Index  to  NASA  Teclmical  Publications 

1957-1961 

NACA  Index  to  NACA  Technical  Publications 

I95O--I956 

Heat  Transfer  Revie'f/s 

Industrial  and  Engineering  Chemistry  (1948-1962) 

Mechanical  Engineering  (196O-I962) 

International  Journal  of  Heat  and  Mass  Transfer  (I96O-I963) 
Including  Heat  Transfer  Bibliography 

American  Rocket  Society  Journal  (1954-1963) 

(Now  AIAA  Journal) 

Including  Technical  Literature  Digest 

Cryogenics  (I962-I965) 

Including  Bibliography  on  Heat  Transfer 


Journal  of  Heat  Transfer  (ASME)  (1959-February  1964) 

Advances  in  Cryogenic  Engineering 
Vol.  1-8  (1955-1962) 

International  Developments  in  Heat  Transfer 

Papers  of  1st  International  Heat  Transfer  Conference  I96I;  ASME 

Jet  Propulsion  Laboratory 

Publications  of  the  Jet  Propulsion  Laboratory  (July  1958-June  1962) 
Applied  Cryogenic  Engineering 

Vance,  R.  W.  and  V.  M.  Duke,  Ed.  John  Wiley  &  Sons  Inc.,  New 
York,  1962 

A  Survey  of  the  Literature  on  Heat  Transfer  from  Solid  Surfaces  to 
Cryogenic  Fluids 

Richards,  R^  J, ,  W.  G.  Steward,  and  R.  B*  Jacobs 
National  Bureau  of  Standards 
Technical  Note  122 
October  I96I 

A  Bibliography  of  Low  Temperature  Engineering  and  Research 
Mendelssohn,  K*,  R.  B.  Scott,  and  L.  Weil 

1944-1960 

Heywood  and  Co.  Ltd.,  London 

Published  as  a  special  issue  of  "Cryogenics" 

A  Brief  Review  of  the  Literature  on  Boiling  Heat  Transfer 
Hawkins,  G.  A. 

University  of  California,  Department  of  Engineering,  Los  Angeles 
June  1950 

Problems  in  Heat  Transfer  During  a  Change  of  State 
A  collection  of  articles 
Kutateladze,  S.  S. 

State  Power  Press,  Moscow,  1953 
(AEC-t rans 1 at i on-AEC-tr-3405 ) 

Current  Contents  (1962-January  1964) 
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Papers  from  9th  Annual  Cryogenic  Engineering  Conference 
Boulder,  Colorado  (August  19^)3) 

Papers  from  ASME-AIChE  Heat  Transfer  Conference 
Boston,  Massachusetts  (August  1963) 

A  total  of  7^1  applicable  references  were  found,  including  replies  to 
letters  sent  to  various  authors  requesting  tabulated  data.  These 
references  are  listed  alphabetically  according  to  author  in  the  Experi¬ 
mental  Data  Bibliography.  It  is  believed  that  the  literature  search  is 
the  most  complete  available  for  boiling  heat  transfer  to  cryogenic  fluids 
from  solid  surfaces. 

Special  attention  was  given  to  articles  in  foreign  languages.  Informal 
oral  translations  were  obtained  from  experienced  scientists  and  engineers 
to  determine  the  content  of  each  of  these  references. 

Of  the  74  references,  62  were  in  English  and  12  were  in  foreign  languages 
as  follows: 

Reference 

Language  No. 

Czechoslovakian  6l 

French  4,  58,  68,  69, 

70,  71,  72,  73 

German  28 ,  29 

Russian  65 

It  is  interesting  to  note  that  although  the  Russians  have  devoted  much 
effort  to  the  development  of  basic  boiling  theories,  they  have  sho^m 
little  apparent  interest  in  obtaining  fimdamental  experimental  data  for 
boiling  cryogenics.  The  only  Russian  papers  found  (Ref.  48  and  65)  were 
concerned  with  obtaining  rather  specialized  boiling  data. 


10 


CODING  SCHEME 


Because  of  the  numerous  geometric  configurations  and  operating  conditions 
for  which  experimental  boiling  data  may  be  obtained,  it  was  considered 
desirable  to  devise  a  coding  scheme  to  allow  rapid  determination  of  the 
boiling  conditions  for  a  given  reference.  Accordingly,  73  of  the  7^ 
references  were  coded  as  given  in  the  Experimental  Data  bibliography. 
Reference  64  was  not  coded  because  it  contains  only  over-all  coefficients 
for  a  condensing- boiling  system. 

The  coding  scheme  developed  is  given  in  Table  1„ 

As  part  of  the  coding  scheme,  each  reference  was  also  rated  with  respect 
to  the  extent  to  which  the  experimental  technique  was  described.  The 
following  rating  scheme  was  used: 

X  :  Little  or  no  discussion  of  experimental  technique 

XX  :  Some  discussion  of  experimental  technique,  but 

insufficient  detail  eonceming  preparation  of 
surface  and/or  instrumentation 

XXX  :  Relatively  complete  details  of  experiments 

This  rating  was  applied  to  the  entire  individual  reference  and  was  placed 
to  the  right  of  the  coding  for  the  boiling  conditions. 
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TABLF  1 


CODING  SCIIENffi 


A.  Boiling  Regime 

1.  Nucleate 

2.  Transition  (or  partial  film) 

3.  Stable  Film 

4.  Peak  Nucleate  (or  burnout  condition) 

B.  Flow  Condition 

1.  Pool 

2.  Forced  Convection — subcooled  (surface  boiling) 

3.  Forced  Convection — saturated  (bulk  boiling) 

4o  Natural  Circulation 


C. 


Ge  ometry 

1.  Above  Horizontal  Plates 

2.  Below  Horizontal  Plates 

3.  Vertical  Plates 

4.  Inside  Vertical  Tubes 

5.  Outside  Vertical  Tubes 

6.  Inside  Horizontal  Tubes 

7.  Outside  Horizontal  Tubes 

8.  Horizontal  Wires 

9.  Vertical  Wires 

10.  Inclined  Plates 

11 „  Miscellaneous  Geometry 
Uc  Circular  Passage 

b.  Annular  Passage  > 

c.  Other  Shaped  Passage 


Pertains  to  items 
4  and  6  only 
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TABUE  1 


(Continued) 

D.  Other  Pertinent  Conditions 

1.  Surface  Effects 

2.  Vortices 

3.  G-field 

4.  Agitation 

5.  Electrical  Field 

E.  Data  Presentation 

1,  Tabulated  Only 

2,  Plotted  Only 

3.  Plotted  and  Tabulated 

4.  Complete  Data  Available  From  Another  Listed  Source 
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An  example  of  the  application  of  the  coding  scheme  as  applied  to  Rel  .  15 
is  as  follows: 

Class,  C.  R.  ,  De  Haan,  J,  R.  ,  Piccone,  M.  ,  and  Cost,  R.  :  ’’Pool 
Boiling  Heat  Transfer  to  a  Cryogenic  Liquid,"  Beechcraffc  Research 
and  Development,  Inc.,  Engineering  Dept.  6l5'i,  WADC  TR  58-5^8, 
AD214256,  1958 

H2:A1-B1-C1,C5,C10-*D1-E5 
H2:A3-B1-C1,C3,C10-D1-E5  XXX 

H2:A4~B1-C1,C3,C10-D1~E3 

From  the  coding,  it  is  immediately  apparent  that  this  work  was  concerned 
with  nucleate,  peak  nucleate,  and  stable  film  pool  boilingo  The  heating 
surface  was  a  flat  plate  in  horizontal,  vertical,  and  inclined  positions^ 
Surface  condition  effects  were  studied,  and  the  data  are  both  plotted 
and  tabulatedo  Relatively  complete  details  of  the  experiments  were 
given. 

Detailed  study  of  the  74  references  in  conjunction  with  the  coding 
scheme  indicated  that,  in  a  number  of  cases,  the  same  experimental  study 
was  reported  in  more  than  one  reference.  The  cases  where  this  occurred 
are  listed  in  Table  2.  It  was  determined  that  the  73  coded  references 
corresponded  to  only  43  separate  experimental  studies. 

The  over-all  scope  of  the  general  conditions  covered  in  the  43  separate 
cryogenic  boiling  studies  is  given  in  Table  3.  It  is  evident  that 
experimental  data  are  completely  lacking  in  a  number  of  cases.  This 
lack  of  data  in  certain  areas  is  discussed  in  detail  in  the  subsequent 
section  on  Recommendations. 
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TABLE  2 


EXPERIMENTAL  STUDIES  REPORTED  IN  MORE  THAN  ONE  REFERENCE 


Reference  Numbers 

Principal  Invest igat or (s) 

1,  56,  62 

0*Hanlon  and  Sherley 

2,  3 

Banchero,  Barker  and  Boll 

5,  6,  7 

Bromley 

8,  9,  10,  11,  12,  13, 

14,  50,  51 

Clark  and  Merte 

15,  16 

Class,  DeHaan,  Piccone  and  Cost 

20,  21 

Dean  and  Thompson 

22,  23 

Drayer  and  Timmerhaus 

24,  25 

Fenster  and  Clark 

CM 

f— I 

CM 

Frederking 

28,  29,  30 

Frederking 

34,  59 

Giaque 

35,  39,  40 

Graham,  Hendricks  and  Hsu 

41,  42 

Hsu  and  Westwater 

44,  66 

Lewis 

46,  47 

Lyon,  Kosky  and  Harmon 

60,  61 

Ruzicka 

67,  74 

Wright  and  Walters 

68,  70 

Weil 

L5 


TABU*  3 


OVER-ALL  SCOPE  OF  CONDITIOXS  COVERED  IN 
EXPERIMENTAL  CRYOGENIC  BOILING  STEDIES 


Conditions 


Nucleate  Boiling: 

Pool 

Subcooled  Forced  Convection 
Saturated  Forced  Convection 

Stable  Film  Boiling: 

Pool 

Subcooled  Forced  Convection 
Saturated  Forced  Convection 

Maximum  Heat  Flux: 

Pool 

Forced  Convection 


Numb 

er  of  Studies 

Hydrogen 

Nitrogen 

Oxygen 

7 

20 

II 

2 

2 

0 

5 

4 

0 

5 

14 

3 

2 

1 

0 

4 

I 

0 

3 

II 

2 

I 

2 

_ _ _ 

0 

In  addition  to  the  alphabetical  listing  by  author  in  the  Experimental 
Data  Bibliography,  a  large  number  of  convenient  cross-reference  bibliog¬ 
raphy  lists  were  prepared.  This  was  facilitated  by  typing  the  references 
on  Keysort  cards,  using  the  coding  scheme  in  punching  the  cards.  These 
listings,  by  reference  number  only,  allow  rapid  determination  of  the 
references  applicable  to  specific  conditions  of  boilingo  The  cross- 
reference  lists  are  included  after  the  Experimental  Data  Bibliography. 
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PHYSICAL  PROPERTIES 


In  connection  with  the  evaluation  of  the  various  theoretical  boiling 
predictions  to  be  presented  in  the  next  section,  an  extensive  compilation 
of  pertinent  thermodynamic  and  transport  properties  of  hydrogen,  nitrogen, 
and  oxygen  was  completed*  Saturation  properties  were  compiled  from  the 
melting  point  to  the  critical  temperature*  Superheated  vapor  properties 
were  tabulated  from  the  normal  boiling  temperature  to  3000  R  and  for 
pressures  from  one  atmosphere  to  the  critical  pressure.  The  melting, 
normal  boiling,  and  critical  conditions  for  the  three  cryogenic  fluids 
are  listed  in  Table  4. 


TABLE  4 

MELTING,  NORMAL  BOILING,  AND  CRITICAL  CONDITIONS  FOR 
HYDROGEN,  NITROGEN,  AND  OXYGEN 


Property 

Hydrogen 

Nitrogen 

Oxygen 

Melting  Temperature,  R 

25.0 

113.8 

98.5 

Normal  Boiling  Temperature,  R 

36.5 

139.2 

162.3 

Critical  Temperature,  R 

59.4 

227.3 

278.6 

Critical  Pressure,  psia 

187.7 

492.9 

736.0 

The  following  saturation  properties  were  compiled: 


Liquid  Density 
Vapor  Density 

Liquid  Specific  Heat  at  Constant  Pressure 
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Liquid  Viscosity 

Liquid  Thermal  Conductivity 

Vapor  Pressure, 

Surface  Tension 

a 

Heat  of  Vaporization 

The  following  superheated  vapor  properties  were 

compiled: 

Vapor  Density 

Pv 

Vapor  Specific  Heat  at  Constant  Pressure 

s 

Vapor  Viscosity 

Vapor  Thermal  Conductivity 

All  properties  are  presented  in  graphical  form  in  Appendix  A.  A  i 
sion  of  the  methods  employed  for  the  compilation  of  each  property 

given  below.  The  references  are  included  in  the  Physical  Property 
References  section. 


DENS  ITT,  AM) 


Hydrogen 

Values  of  the  density  of  saturated  liquid  and  saturated  vapor  from  35  P- 
to  the  critical  temperature  were  taken  from  the  tabulated  values  of 
NBS  TN-I30  (Ref.  101).  Values  down  to  25  R  were  extrapolated  from  a 
plot  of  the  above  values.  The  superheated  vapor  densities  were  calcu¬ 
lated  from  various  equations  of  state,  as  summarized  below. 
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Temperature  Range , 
_ R _ 


Reference 


36-180 

180-540 

540-2500 

2500-3000 


102 

101 

104 

perfect  gas  relations 


Nitrogen 

All  saturation  properties,  both  liquid  and  vapor,  vere  taken  from 
Johnson  (Ref.  104).  The  superheated  vapor  values  to  1000  R  were  calcu¬ 
lated  from  the  equation  of  state  recommended  by  NBS  (Ref.  105)«  Above 
this  temperature,  the  perfect  gas  relationships  were  assumed  to  apply. 


Oxygen 

Values  of  saturated  liquid  density  from  110  to  275  R  were  taken  from 
Johnson  (Ref.  104).  The  lower  temperature  values  were  extrapolated  from 
those  values.  The  values  of  saturated  vapor  density  were  taken  from 
the  Battelle  Handbook  (Ref.  106)«  Superheated  vapor  values  to  1000  R 
were  calculated  from  the  equation  of  state  recommended  by  NBS  (Ref.  10?) . 
From  1000  to  3000  R  the  perfect  gas  relations  were  assumed  to  be 
applicable. 


21 


SPECIFIC  HEAT  AT  CONSTANT  PRESSUEE,  C.  .VNI)  C,. 

Ju  V 

Hydrogen 

The  values  of  saturated  liquid  specific  heat  at  constant  pressui'o  were 
taken  from  Johnson  (Ref,  104)  j  except  for  a  small  extrapolation  do^vTl  to 
25  R»  The  superheated  vapor  values  to  2^00  R  were  taken  from  Ref,  102, 
Above  this  temperature,  values  wort'  taken  from  Hi'f,  103. 


Nitrogen 

The  saturated  liquid  values  from  115  to  210  R  were  taken  from  Johnson 
(Ref.  104).  Above  this  temperature,  the  above  data  were  graphically 
extrapolated.  The  superheated  vapor  values  to  1000  R  were  calculated 
by  the  NBS  recommended  equations  (Ref.  105)*  Above  this  temperature, 
the  values  were  taken  from  Ref.  113.^^ 


^gen 


Specific  heat  values  of  the  saturated  liquid  from  105  to  I60  R  were  taken 
from  Johnson  (Ref,  104).  Above  I60  R,  no  values  could  be  obtained. 
However,  values  of  ,  the  heat  capacity  of  the  liquid  along  the 
saturation  line,  were  available  (Ref,  108).  The  values  of  C  were 

li 

converted  to  by  the  method  of  Hougen  and  Watson  (Ref,  IO9).  The 
specific  heat  of  the  saturated  vapor  was  taken  from  Ref.  107  up  to 
1000  R.  Above  this  temperature,  the  values  were  taken  from  Ref.  II3. 
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VISCOSITY,  IJL^  AND 


Hydrogen 

Values  of  saturated  liquid  viscosity  from  30  to  35  R  were  taken  from 
Johnson  (Ref.  104).  The  remainder  of  the  saturated  liquid  values  and 
the  superheated  vapor  values  were  obtained  from  the  residual  method  first 
suggested  by  Abas-Zade  (Ref.  IIO).  This  method  states  that  the  viscosity 
is  composed  of  an  ideal  low-pressure  contribution,  which  is  a  function 
of  temperature  only  and  a  residual  value  due  to  pressure  m  The 

residual  is  a  function  of  pressure  alone.  The  values  vere  taken  from 
Starrett  (Ref.  Ill)  and  the  (jLir-/i*)  values  from  the  correlation  of 
Thodos  (Ref.  112). 


Nitrogen 

Values  of  saturated  liquid  viscosity  from  115  to  210  R  were  taken  from 
Johnson  (Hef.  104).  From  210R  to  the  critical  point,  the  values  were 
extrapolated.  The  superheated  vapor  values  were  obtained  by  the 
residtxal  method.  The  values  were  taken  from  NBS  circular  564 
(Ref.  113)  and  (/X-pi*)  values  from  Thodos  (Ref.  114). 


Oxygen 

The  values  of  saturated  liquid  viscosity  were  taken  from  Johnson  (Ref.  104). 
The  superheated  vapor  values  were  calculated  by  the  residual  method  with 
^  values  from  Ref.  113  and  (/li-^*)  from  Ref.  114. 
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THERMAL  CONDUCTIVITY,  k.  AKD  Ic. 

Jj  \ 


Hydrogen 

All  values  of  thermal  conductivity  vere  taken  from  the  reduced  state  plot 
of  Thodos  (Rof.  115),  where  is  plotted  against  with  as  a 
parameter. 


Nitrogen 

Values  of  saturated  liquid  thermal  conductivity  from  123  to  160  11  were 
taken  from  Jolmson  (lief.  104).  The  remainder  of  tlie  liquid  values  were 
obtained  from  Thodos  (Uc.''".  113)  ^nd  smoothed  to  fit  the  tabulated  data 
of  Johnson  (Ref,  104) ^  The  superheated  vapor  thermal  conductivity  was 
calculated  by  the  residual  method,  witli  k'-  taken  from  Ref,  113  ai^d 
(k-k*)  from  Ref o  ll6. 


Oxygen 

The  values  of  liquid  thermal  conductivity  from  IpO  R  to  the  critical  point 
were  taken  from  Johnson  (jlef.  104).  The  values  from  100  to  123  R  were 
extrapolated.  Superheated  vapor  thermal  conductivity  was  obtained  by 
the  residual  method,  with  k*  from  Ref.  lip  and  the  residual  from  Ref.  Il6, 
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VAPOR  PRESSURE,  P^ 

Hydrogen,  Nitrogen,  Oxygen 

The  vapor  pressure  of  hydrogen,  nitrogen,  and  oxygen  over  the  entire 
temperature  range  was  obtained  from  Johnson  (Ref,  104). 

SURFACE  TENSION,  O 

Hydrogen 

The  surface  tension  of  hydrogen  was  obtained  from  NBS  7616  (Ref*  11?) • 


Nitrogen 

The  surface  tension  of  nitrogen  was  obtained  by  the  equation  of 
Guggenheim  (Ref.  118)  where: 

a  =  cr  (1  -  t/t  )^V9 

O  ^  '  C 


and  (T^  is  a  constant  for  a  given  fluid*  The  value  of  CT^  was  determined 
at  the  normal  boiling  point  value  given  by  Johnson  (Ref.  104). 


Oxygen 

The  surface  tension  of  oxygen  was  obtained  in  the  same  manner  as  that 
for  nitrogen* 
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HEAT  OF  VAPORIZATION,  X 


Hydrogen 

The  heat  of  vaporization  of  hydrogen  from  55  R  to  the  critical  was  taken 
from  the  tabulated  values  in  KBS  TNI30  (Ref.  101 ).  Values  down  to  25  R 
were  obtained  from  Roder,  et  al.  (Ref.  119)* 


Nitrogen 

The  heat  of  vaporization  of  nitrogen  over  the  entire  saturation  curve 
was  obtained  from  Johnson  (Ref.  104). 


Oxygen 

The  heat  of  vaporization  of  oxygen  was  obtained  from  the  Watson  equation 

(Ref.  120): 


where  is  the  heat  of  vaporization  at  any  reduced  temperature,  T 
The  reference  value  of  X^^  'Vvus  taken  at  the  atmospheric  boiling  ^ 
point  from  Ref.  120. 
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BOILING  THEORY 


Prior  to  1950,  only  a  few  empirical  correlations  of  boiling  data  were 
available o  Since  then,  however,  a  considerable  number  of  theoretical 
and  semi theoretical  equations  have  been  developed  to  predict  boiling 
phenomena  in  both  the  nucleate  and  stable  film  boiling  region,  as  well 
as  at  the  burnout  point.  This  has  been  particularly  true  for  pool 
boilingo  This  is  illustrated  in  Fig.  2  where  a  complete  boiling  curve 
for  a  specific  case  of  hydrogen  is  constructed  solely  on  the  basis  of 
theories  to  be  discussed  below.  Some  notable  theoretical  progress  has 
also  been  made  for  subcooled  forced  convection  (surface)  boilings  How¬ 
ever,  for  saturated  forced  convection  (bulk)  boiling  and  the  accompany¬ 
ing  complex  two-phase  flow  phenomena,  only  a  start  has  been  made  and 
little  is  known. 

This  section  presents  a  review  of  all  of  the  available  theoretical 
approaches  to  boiling  phenomena.  In  all  cases  the  theories  are  general 
and  are  not  restricted  to  cryogenic  fluids o  However,  applications  of 
the  equations  in  this  report  are  made  only  to  the  cryogenic  fluids  of 
interest.  In  general,  empirical  expressions  are  not  considered. 
Included  are  equations  covering  the  following  areas  of  boiling: 

Pool  Nucleate  Boiling 

Pool  Nucleate  Boiling  Maximum  Heat  Flux 
Pool  Film  Boiling  Minimum  Heat  Flux 
Pool  Stable  Film  Boiling 

Forced  Convection  Subcooled  Nucleate  Boiling 

Forced  Convection  Subcooled  Nucleate  Boiling  Maximum  Heat  Flux 
Forced  Convection  Saturated  Boiling 
Forced  Convection  Film  Boiling 
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Wax. 

Nucleate 
(Eq.  28) 


FIGURE  2  TYPICAL  THEORETICAL  POOL  BOILING  CURVE,  HYDROGEN  -  12  psia  -  HORIZONTAL  PLATE  PACING  UP 


A  fundamental  discussion  of  boiling  is  not  included  in  this  report,  but 
may  be  found  in  a  number  of  excellent  references  including  McAdams 
(Ref.  202),  Jakob  (Ref„  203),  Grober,  Erk  and  Grigull  (Ref,  204), 

Rohsenov  (Refc  205),  and  Westvater  (Ref,  206).  For  the  most  part,  com¬ 
parison  of  theory  with  experimental  data  is  discussed  in  the  next  section 
of  this  report. 


POOL  NUCLEATE  BOILING 

To  initiate  pool  boiling,  the  heat  transfer  surface  temperature  must  be 
somewhat  higher  than  the  saturation  temperature  of  the  fluid  at  the 
existing  pressure.  If  the  surface  temperature  is  less  than  the  value 
required  to  initiate  boiling,  then  heat  transfer  will  occur  by  natural 
convection  only.  This  process  is  governed  by  equations  similar  to  the 
following  for  natural  convection  above  a  horizontal  flat  plate,  as  given 
by  McAdams  (Refo  202): 


Thus,  the  heat  flux  is  proportional  to  the  temperature  difference  raised 
to  an  exponent  somewhat  greater  than  one.  Nonboiling  relations  for 
other  geometrical  configurations  are  also  given  in  Ref,  202  and  are  not 
discussed  here. 

When  a  sufficient  superheat  exists  to  cause  nucleate  boiling  then,  as 
shown  in  Fig.  2,  the  heat  flux  becomes  a  stronger  function  of  the  tem¬ 
perature  difference.  With  increasing  A  T,  the  bubble  population  grows 
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until  it  reaches  a  maximum  density  corresponding  to  the  maximum  pool 
nucleate  boiling  heat  flux.  For  the  AT  region  between  the  initiation 
of  boiling  and  this  maximum,  a  number  of  theories  have  been  proposed 
for  saturated  pool  nucleate  boiling.  These  are  summarized  in  Table  5 
as  Eq.  2  through  10.  In  general,  these  equations  can  be  arranged  into 
the  following  general  boiling  Nusselt  number  form  suggested  by  Rohsenow 
(Ref,  301)  and  discussed  by  Zuber  and  Fried  (Ref^  302); 

=  C,  Re  ™  Pr^  (ll) 

D  1  D 


or,  tetter,  in  an  equivalent  boiling  Stanton  number  form  to  separate 
the  heat  flux,  q/A,  from  the  superheat,  AT: 


Reg 

^2 

Stg"  Pr^ 

(12) 

where : 

A  T 

=  T  -  T 
w  3 

(13) 

NUr 

hL* 

(14) 

L*G* 

(15) 

Pr 

II 

0 

tr' 

(16) 

CrG*  (At)CrG* 

(17) 

In  effect,  the  boiling  Stanton  number  is  a  ratio  of  the  latent  heat 
effect  to  the  sensible  heat  effect.  The  characteristic  length,  L*,  in 
the  Reynolds  number  is,  in  most  cases,  based  on  an  initial  bubble  radius. 
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the  bubble  diameter  at  departure  such  as  in  Eq.  2,  3,  and  7,  a  bubble 
diameter  referred  to  some  time  constant  such  as  in  Eq.  4,  or  a  vapor 
front  advance  distance  such  as  in  Eq.  8.  The  characteristic  superficial 
mass  velocity,  G*,  in  the  Reynolds  number  is,  in  most  cases,  based  on 
the  liquid  assuming  either  Gj^=Gy  as  in  Eq.  2,  or  as  in  the  other 

equations  of  Table  5.  Table  6  compares,  in  Stanton  number  form,  the 
equations  listed  in  Table  5.  In  developing  Table  6,  the  folloving 
assumptions  were  made  to  simplify  the  relations: 


,  (A  T)  X 

Ap  -  j  '' 

S 

(18) 

(pl  -  'V>  “ 

(19) 

In  all  cases,  G*  (q/A).  Thus,  the  Reynolds  number  Re^  “  (q/A)  and 
l/Stg  “  (At)®,  thus  separating  out  the  effects  of  q/A  and  A  T  as  noted 
above. 

In  general,  the  theories  predict: 

q/A  “(At)"  2  ^  n  ^  3-l/3 

and  q/A  «  (^^L^  ^ 

with  the  exponent,  -b,  ranging  from  -5.1  to  -19/24.  The  effect  of  liq¬ 
uid  viscosity  cancels  out  in  the  equations  of  Levy  (Ref.  306),  Michenko 
(Ref.  307),  Chang  and  Snyder  (Ref.  309)  and  Nishikawa  and  Yamagata 
(Ref.  310).  The  predicted  variation  of  the  heat  flux  with  the  surface 
tension  is  given  by: 

q/A  “  (a)'^ 

with  c  ranging  from  -I0/6  to  -I/2. 
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The  pool  nucleate  boiling  equations  listed  in  Table  5  vere  programmed 
for  the  IBM  7094  to  facilitate  calculations  for  the  cryogenic  fluids  of 
interest . 

Values  of  the  required  physical  properties  for  N^,  and  O^,  given  in 

Appendix  A,  were  included  in  the  computer  program  to  permit  direct  appli¬ 
cation  of  the  equations  to  these  fluids.  Results  of  calculations  for 
a  pressure  of  about  1  atmosphere  are  given  in  Fig.  4,  5,  and  6  for  hydro¬ 
gen,  nitrogen,  and  oxygen,  respectively.  A  value  of  =  0.015  fop 
the  Rohsenow  equation  was  used.  Calculations  were  also  carried  out  for 
a  wide  range  of  pressure.  Results  from  the  Forster  and  Zuber  equation 
are  shown  in  Fig.  7,  8,  and  9.  This  equation  appeared  to  give  results 
representative  of  the  better  nucleate  boiling  equations.  The  increase 
in  heat  flux  with  increasing  pressure  for  a  given  A  T  is  typical  of  the 
trends  of  experimental  data  as  is  discussed  in  the  following  section. 

The  curves  in  Fig.  7,  8,  and  9  are  terminated  at  approximately  the  pool 
nucleate  boiling  maximum  heat  fluxes  as  determined  from  the  limited 
experimental  data  of  Class,  et  al.  (Ref.  15  and  l6)  for  hydrogen  and 
Lyon,  et  al.  (Ref.  46)  for  nitrogen  and  oxygen  for  boiling  above  hori- 
zontal  plates „ 

In  general,  the  greatest  shortcoming  of  the  above  theories  for  pool 
nucleate  boiling  is  their  failure  to  theoretically  account  for  variables 
which  are  not  easily  contr oiled „  These  include  amounts  of  dissolved  and 
adsorbed  gases,  and  impurities,  contamination,  and  nature  of  the  heating 
surface o  Also  of  importance  are  the  size,  geometry  and  orientation  of 
the  heating  surface.  This  is  discussed  more  fully  in  the  next  section 
where  experimental  data  are  presented.  These  effects  can  influence  both 
the  slope  and  position  of  the  boiling  curve.  Consequently,  the  equations 
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SLMMAEY  OF  POOL  NUCLEATE  BOILING  EQUATIONS 


Reference 


Rohsenow 

(301) 


Kutateladze 

(303) 


Forster  and 
Zuber  (304) 


Equation 


Cl  (At) 


=  c 


fiaAi 


Sf  I  ^  I  g  (Pl  -  Py) 


gc  cy 


]"} 


0.35 


1.7 


Typically  C  „  varies  from  0.0027  to  0.015  depend 

S  X 

ing  upon  fluid-heater  surface  combination. 


n  1/2 


L  ®  (^L  - 


=  7.0  X  10 


-4 


r Pl  (^y^) 

1  ^  j_  g  (Pl  ”  Py^  J 

(ClVV'-'' 


1/2  0.7 

} 


0.70 


TABLE  5 
(Continued) 


Reference 

Equation 

Equation 

No. 

Forster  and 
Greif  (305) 

1/0 

q/A  =  0.0012  (k^  T^/j  X  a  ^  )  * 

(L L  h/-’  • 

5 

Levy  (306) 

Jk  c  p^(At)^ 

i/a  -  a  (p,  -  ,v)  (Pv-  « 

The  function  (p^  X)  is  given  in  Fig.  3 

6 

Michenko 

(307) 

_  ^  1/2 

^  —4 

■  p  7  Y  1  n 
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TABLE  5 
(Continued) 


Reference 

Equation 

Equation 

No. 

Labountzov 

(308) 

'“’l  ^  ■'s  ,  r  <^1  Pl  ''  "s 

8 

kj^  (p^  X)^  J  ^  {p^xf  J  \  h  ) 

V - ^ - / 

Re 

—2 

For  Re  >  10  ,  m  =:  0c65,  C,  =  0ol25 

1 

For  Ee  <  10“^,  m  =  0»50,  =  0.0625 

Chang  and 

Snyder 

(309) 

h  =  4  X  10  * 

(Pl 

1  0.4 

-  Py) 

9 

Nishikawa 

and 

Yamagata 

(310) 

fi=  8.0- 

‘‘l 

m  =  274.32 

/•' 

< 

f 

Pl  Pl'  s 

B  kL  a  X  Py  go 

t“^  B  =  6.742  I 

1/2  1  2/3 

(q/A)  >  L  • 

-/ 

5tu/hr 

10 

Heat  Flux,  q/A,  Btu/hr-ft 
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PREDICTED  POOL  NUCLEATE  BOILING  HEAT  FLUX  FOR 
NITROGEN  ACCORDING  TO  FORSTER  AND  ZUBER  EQUATION 
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PREDICTED  POOL  NUCLEATE  BOILING  HEAT  FLUX  FOR  OXYGEN 
ACCORDING  TO  FORSTER  AND  ZUBER  EQUATION 


are  not  general  and  cannot  be  expected  to  predict  the  heat  transfer  con¬ 
dition  accurately  for  any  heater-fluid  combination.  Some  recent  work 
on  improved  theories  which  accounts  for  the  nature  of  the  heating  sur¬ 
face  and  its  relationship  to  active  sites  has  been  presented  by  Kurihara 
and  Myers  (Ref.  311)5  Tien  (Ref.  312),  and  Lienhard  (Ref,  313).  For 
example,  the  latter  equation  is  of  the  form: 


q/A 


r  org  g  (p  n  )  1  1/2  5/4 

V  -I 


1/3 


The  number  of  active  sites,  n,  is  a  function  of  the  nature  of  the  sur¬ 
face  and,  in  general,  increases  as  the  temperature  of  the  surface 
increases. 


Furthermore,  theoretical  development  of  pool  nucleate  boiling  for  cryo¬ 
genic  fluids  will  also  be  limited  unless  detailed  boiling  mechanism 
studies,  specifically  for  cryogenic  fluids,  are  conducted.  An  initial 
start  in  this  important  area  has  been  made  by  McFadden  and  Grassman 
(Ref.  314).  Their  study  concerned  bubble  frequency  and  diameter  at  a 
given  nucleation  site.  Experimental  data  were  obtained  with  liquid 
nitrogen  and  were  compared  with  the  correlation  of  Zuber  (Ref.  323): 

f  D  =  0.59  - ^ 
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This  correlation  had  been  shown  to  work  reasonably  well  for  water,  car¬ 
bon  tetrachloride,  and  methyl  alcohol.  The  nitrogen  data,  however,  did 
not  fit  the  Zuber  correlation.  Accordingly,  it  was  shown  by  dimensional 
analysis  that: 


2 

f  D^P 


1/2 


D) 


1/2 


g  (^- 


(21) 


and  a  plot  of  all  data,  including  nitrogen,  was  best  represented'by; 
.  0.56  I  ‘  'VA 


(22) 


Far  from  the  critical  temperature  (fi^-  /Q^)  cw  and  accordingly: 

f  =  0.56  (g)l/2  (23) 


Equation  23  correlated  all  available  results  more  closely  than  Eq.  20; 
however,  the  authors  felt  that  more  experimental  work  is  needed  to 
extend  the  diameter  range  and  a  wider  variety  of  fluids  should  be  studied. 
Relations  such  as  Eq.  23  will  be  useful  in  developing  more  exact  nucleate 
^®ill^g  theories  for  cryogenics  based  on  active  sites. 


POOL  NUCLEATE  BOILING  MAXIMUM  HEAT  FLUX 

The  regime  of  pool  nucleate  boiling  is  limited  to  a  maximum  heat  flux 
and  a  corresponding  maximum  At.  As  with  other  fluids,  the  maximum 
heat  flux  and  maximum  At  for  cryogenic  fluids  both  depend  on  the  par¬ 
ticular  fluid,  the  pressure  and,  to  some  extent,  the  geometry  and 
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orientation  of  the  surface „  The  maximum  heat  flux  is  essentially 
independent  of  the  nature  of  the  surface  (except  if  it  is  greased)  but 
the  maximum  is  highly  dependent  on  this  factor. 

A  number  of  theories  have  been  proposed  recently  to  estimate  the  pool 
nucleate  boiling  maximum  heat  flux.  Some  postulate  that  this  maximum 
occurs  at  the  point  where  the  vapor  formation  rate  is  of  such  magnitude 
that  hydrodynamic  instabilities  begin  to  impede  the  counterflow  of  the 
liquid  toward  the  heating  surface.  All  the  theories  proposed  can  be 
rearranged  easily  in  terms  of  the  superficial  vapor  disengaging  velocity: 


(q/A) 


max 


(24) 


*It  should  be  noted  that  this  is  not  a  dimensionless  parameter^ 


where  includes  buoyancy,  acceleration  and/or  other  physical  property 

effects.  Expressions  for  of  the  theories  reviewed  are  summarized 

in  Table  7.  The  importance  of  certain  liquid  properties  such  as  , 

^  L  ^  major  point  of  contention  among  the 

theories . 


Calculations  were  carried  out  for  the  pool  nucleate  boiling  maximum  heat 
flux  as  a  function  of  pressure  using  the  theories  summarized  in  Table  7. 
The  results  are  plotted  in  Fig,  12,  13,  and  14  for  N^,  and  0^, 
respectively.  In  general,  considerable  disagreement  between  the  various 
theories  exists  over  the  entire  pressure  range.  As  is  typical  with  all 
fluids,  a  maximum  pool  nucleate  boiling  maximum  heat  flux  is  predicted 
at  a  pressure  about  one-third  the  critical  pressure.  Comparison  with 
experimental  data  is  discussed  in  the  next  section. 


•  iiT 


TABLE  7 


SUMMAEY  OF  POOL  NUCLEATE  BOILINTi 
HEAT  FLL^'  TUEOBIES 


(q/A) 


r — (consistent  units  except  as  noted) 
A  ^  X 
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TABLE  7 
(Continued) 
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COKPAKISON  OP  THEORIES  FOR  POOL  NUCLEATE  BOILING  MAXIMUM  HEAT 
FLUX  FOR  HYDROGEN 
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PIG.  13  COMPARISON  OP  THEORIES  POR  POOL  NUCLEATE  BOILING  MAXIMUM  HEAT  FLUX  FOR  NITROGEN 
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FIG.  14-  COMPARISON  OF  ?H30?J3S  FOR  POOL  NUCLEAT3  BOILING  M-AXIMUM 
HEAT  FLUX  FOR  OXYGEN 


POOL  FILM  BOILING  MINIMUM  HEAT  FLUX 


Just  as  the  nucleate  boiling  regime  is  limited  to  a  maximum  heat  flux, 
the  film  boiling  regime  is  limited  to  a  minimum  heat  flux  as  shovn  in 
Fig.  2„  Available  theories  for  the  pool  film  boiling  minimum  heat  flux 
are  summarized  in  Table  8  in  terras  of  the  same  parameter  used  for  the 
maximum  heat  flux,  namely: 

(35) 

All  the  theories  are  based  on  hydrodynamic  instability  of  the  liquid- 
vapor  boundary  and  strictly  apply  only  to  boiling  above  a  horizontal  flat 
surface.  In  all  cases,  the  physical  properties  are  evaluated  at  the 
saturation  temperature,  except  (p^)  which  is  evaluated  at  the  average 
vapor  film  temperature.  ^ 


Calculations  of  the  minimum  heat  flux  were  carried  out  using  the  equations 
in  Table  8  for  each  of  the  three  cryogenic  fluids.  The  results  are  plot¬ 
ted  in  Fig,  15,  16,  and  I7  for  H2,  N^,  and  0^,  respectively.  By  using 
the  parameter  1  ordinate  rather  than  just 

possible  to  plot  the  theories  in  terms  of  pressure  only. 
At  relatively  low  values  of  At,  these  two  vapor  densities  approach  each 

other  and  the  ordinate  approaches  the  value  of  (q/a)  At  larirer  A  T 

'  min  ^ 

values,  the  ordinate  is  somewhat  higher  than  Experimental  data 

for  the  pool  film  boiling  minimum  heat  flux  of  cryogenics  are  very  meager, 
as  discussed  in  the  next  main  section. 


TABLE  8 


SUMMARY  OF  POOL  FIIM  BOILING  MINIMUM 
HEAT  FLUX  THEORIES 
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=  4>2  (consistent  \mits) 
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(Pl  ^ 

1/4 
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36 

Zuber  and 
Tribus  (318) 
(alternate 
approach) 

T  1/4  K„  =  0.109  to  0.144 

a  gg^  ^  2 

Pl  “  Py 

37 

Zuber  (324) 
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Berenson 
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FIG.  15  POOL  FILM  BOILING  MINIMUM  HEAT  FLUX  FOR  HYDROGEN 
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POOL  FILM  BOILING  MINIMUM  HEAT  FLUX  FOR  OXYGEN 


POOL  STABLE  FILM  BOILING 


The  pool  stable  film  boiling  regime  is  important  for  cryogenic  fluids 
because  it  can  be  achieved  at  surface  temperatures  within  the  practical 
limits  of  a  number  of  materials.  In  fact,  it  is  possible  to  attain 
higher  heat  fluxes  than  the  pool  nucleate  boiling  maximum  heat  flux  by 
operating  under  conditions  of  film  boiling.  This  regime  is  much  more 
amenable  to  theoretical  development  than  pool  nucleate  boiling.  Since 
the  original  theoretical  derivations  by  L.  A,  Bromley  in  1948  (Ref.  5  and 
6),  based  on  the  similarity  to  film  condensation,  a  great  deal  of  effort 
has  been  expended  on  improving  the  theory  and  extending  it  to  different 
boiling  surface  geometries.  The  theories  are  best  summarized  with  respect 
to  the  geometric  configuration.  Corrections  for  the  effect  of  thermal 
radiation  through  the  vapor  film  at  high  temperatures  may  be  important 
and  are  discussed  at  the  end  of  this  section.  It  is  important  to  note 
that  nitrogen  was  one  of  the  fluids  used  by  Bromley  initially  to  verify 
his  theory. 


Long  Horizontal  Cylinders 


Bromley’s  original  equation,  assuming  laminar  flow,  was: 


h 

c 


Py  iPi;-Py)  g 

^  K  -  V  ^ 


(40) 


where;  The  constant  C^  was  found  to  be  bracketed  theoretically  from 
0.512  (stagnant  liquid  surrounding  the  tube)  to  0.724  (liquid 
moves  completely  freely  with  the  vapor).  An  average  experi¬ 
mental  value  of  C^  =  0,62  was  determined  by  Bromley  from 
experimental  data. 


60 


and : 


X'  =  X  t  0.5  (T^  -  T^) 


(41) 


All  vapor  properties  except  in  (p^  -  p^)  terras  are  evaluated  at: 


T 

a  ve 


T  +  T 
V  s 


2 


Bromley  later  (Ref.  326)  modified  Eq.  41  to  take  into  account  more 
properly  the  superheating  of  the  vapor  film: 


=  X 


1  +  0.4 


V 


(T  - 
^  w 


T  ) 
s 


(42) 


Rohsenow  (Ref.  32?)  further  modified  Eq .  42  to  account  for  crossflow 
within  the  film: 


1  +  0.675 


C,,  (t  -  T  ) 
V  ^  w  s^ 


X 


1  -h  0,34 


X 


(43 


\ 


More  recently,  McFadden  and  Grosh  (Ref.  328)  determined  the  effects  of 
variable  density  and  specific  heat  for  the  vapor  film  by  numerically 
solving  the  boundary  layer  equations  for  stable,  free  convection  film 
boiling.  They  considered  a  vertical  surface  as  we'  1  as  a  horizontal 
cylinder.  They  found,  in  general,  that  compared  to  the  assumption  of 
constant  properties,  a  variable  specific  heat  tends  to  increase  the 
heat  transfer  coefficient  while  a  variable  density  tends  to  thicken  the 
vapor  layer  and  thereby  decrease  the  heat  transfer  coefficient.  The 
latter  effect  was  predominant,  but  in  general,  only  at  pressures  apx“r?. e/ 
ing  the  critical  (about  >  0,9)  the  influence  of  variable  properties 


6l 


important.  Frederking  (Ref.  329)  also  obtained  a  numerical  boundary 
layer  solution  for  a  vertical  plate  and  a  horizontal  cylinder.  His 
method  accounts  for  the  influence  of  momentum  transport  at  the  vapor- 
liquid  interface  and  covers  a  wide  range  of  vapor  superheat. 


An  alternate  form  of  Eq.  40  in  terms  of  dimensionless  groups  as  dis¬ 
cussed  by  Frederking  (Ref.  330)  is: 


Nuj^  =  0.62 


1/4 


where : 


h  D 

Nusselt  number  =  - 


Ra^^  =  Rayleigh  number  = 


.  (44) 

(45) 

(46) 


In  general,  Eq.  40  and  44  are  accurate  only  for  a  restricted  range  of 
cylinder  diameter.  For  very  small  diameters,  the  vapor  film  thickness 
may  approach  the  diameter  (this  is  not  allowed  for  in  the  theory)  and 
predicted  values  of  h^  may  be  somewhat  lower.  For  very  large  diameters, 
an  h^  of  zero  is  predicted  which  is  not  consistent  with  experimental 
findings.  Modifications  of  Eq.  40  for  the  effect  of  diameter  over  a 
wide  range  are  discussed  by  Banchero,  et  al.  (Ref.  3)  and 
Frederking  (Ref.  27)  in  the  light  of  experimental  data  for  horizontal 
wires  and  large  cylinders.  Banchero  in  particular,  on  the  basis  of 
extensive  data  for  oxygen,  recommends  that  the  constant  in  Eq,  40  be 
replaced  by  the  expression: 


a  (^  +  c)  D 


1/4 


(47) 


vhere  a  and  c  must  be  determined  from  experimental  data. 
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Horizontal  Plates  (Facing  Up) 


Based  on  a  wave  theory  for  natural  convection,  Chang  (Ref.  33l)  developed 
the  following  equation  which  is  rearranged  in  terms  of  the  Rayleigh  number 

Nuj^  =  0.278 

c 

where:  L  is  substituted  for  D  in  Eq.  45  and  46. 

c 

It  should  be  noted  that  the  characteristic  length,  L  ,  in  Eq.  4?  cancels 

c 

out . 


Ra 


(T  -  T  ) 
V  '  w  s 


1/3 


(47) 


More  recently,  Berenson  (Ref.  332)  used  the  Taylor  instability  concept 
to  derive  the  following  equation  in  terms  of  the  Rayleigh  number  for 
the  region  of  the  minimum: 


Nub 


0.425 


where : 


B 


is  substituted  for  D  in 
Eq.  45  and  46. 


(48) 


(49) 


Vertical  Surfaces 


Bromley  (Ref.  5  and  6)  also  derived  an  expression  for  film  boiling  from 
a  vertical  surface  assuming  laminar  flow.  In  terras  of  the  average 
Nusselt  number  over  a  vertical  length,  L,  and  the  Rayleigh  number,  his 
expression  was:  ^ 


Nu^  =  C, 
L  4 


Ra 


X’ 

L  (T  -  T  ) 
'  V  ^  V  s 


1/4 


(50) 


where:  L  is  substituted  for  D  in  Eq.  45  and  46. 
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The  constant  in  Eq .  50  vas  theoretically  determined  by  IIsu  and  West- 
water  (Ref.  42)  to  be  bracketed  between  O.667  and  0.943.  Bonilla  (Ref,  333) 
suggested  an  average  value  of  0,80. 

Ellion  (Ref.  334)  derived  an  equation  similar  to  Eq.  50  for  the  case  of 
uniform  heat  flux.  He  obtained  a  constant,  of  O.7I6.  More  recently, 

a  number  of  laminar  boundary  layer  solutions  for  vertical  surfaces  have 
been  developed,  including  those  of  McFadden  and  Grosh  (Ref.  328),  Cess 
(Ref.  335),  Sparrow  and  Cess  (Ref.  336),  Koh  (Ref.  337),  Koh  and  Nil son 
(Ref.  338),  and  Frederking  (Ref.  329).  The  aim  of  these  more  refined 
analyses  has  been  to  account  for  inertia  forces,  convective  energy  trans¬ 
port,  interfacial  velocity,  variable  physical  properties  and  interacting 
radiation  and  convection  effects. 


Koh ' s  exact  laminar  boundary  layer  analysis  for  constant  properties,  in 
particular,  results  in  a  lucid  description  of  the  effects  of  shear  stress 
and  vapor  velocity  at  the  vapor-liquid  interface.  He  shows  that  the 
constant  C^  in  Eq.  50  of  Bromley  is  actually  a  function  of: 


-,1/2 


and 


Cy  (At) 


For  values  of 


'V 


1/2 


>  0,001  and  values  of  C^(A  t)/X  (Pr)^, 


Koh  predicts  values  of  h^  considerably  higher  than  those  predicted  by 


Eq.  50. 


In  many  practical  situations,  the  above  laminar  flow  analyses  for  verti¬ 
cal  surfaces  are  not  valid.  Hsu  and  Westwater  (Ref.  338)  observed  heat 
fluxes  100  to  300  percent  greater  than  predicted  by  Eq.  50.  The  increase 
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was  caused  by  the  development  of  turbulent  flow  above  a  critical  height 
corresponding  to  a  critical  vapor  film  thickness  and  film  thickness 
Reynolds  number.  The  critical  Reynolds  number  for  transition  from  vis¬ 
cous  to  turbulent  flow  at  a  height  is  given  by: 


y*  u^*  (p^) 


Re^  = 


100 


where : 


^  ^  Re* 


K),  r  '’l  -  ('v)  1 

f  L  f  J 


1/3 


(52) 


and : 


Py  Re*  X'  y* 

2  k^  (At) 


(53) 


For  Re*  >  100,  they  developed  an  equation  for  h^  which  accounted  for 
the  turbulent  flow  effects.  Bankoff  in  a  written  discussion  to  Ref.  338 
presented  the  following  simpler  alternate  turbulent  flow  equation  for 
the  average  heat  transfer  coefficient: 


h 


-  0.29 


,  V  2/5 

- 

Cy  (At) 

V  '’v  j  1 

— X — J 

ll/5 


V 


(54) 


The  constant  0.29  in  Eq. 
ing  on  the  nature  of  the 


54  is  actually  subject  to  some  variation  depend- 
fluid. 
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Radiation  Effects 


In  his  original  derivation  for  a  horizontal  cylinder,  Bromley  (Ref.  5  and 
6)  also  included  an  expression  for  the  transfer  of  heat  from  the  cylinder 
vail,  through  the  vapor,  to  the  liquid  by  thermal  radiation.  He  noted 
that  the  coefficient  for  convection,  h^,  vas  dependent  on  the  effective 
coefficient  for  radiation,  h^.  For  h^  <  h^,  he  developed  the  following 
approximate  relations: 


total 


he  +  (3/4)h^ 


where : 

o'  (t  ^  -  T 


(55) 


(56) 


In  the  case  of  boiling  on  a  0.350-inch  horizontal  carbon  tube  at 
atmospheric  pressure,  he  found  that  the  radiation  contribution  became 
important  above  a  wall  temperature  of  about  500  F.  More  detailed  evalu¬ 
ations  of  the  effect  of  thermal  radiation  on  film  boiling  have  been  pre¬ 
sented  by  Chang  (Ref.  331 )  and  Koh  and  Nils  on  (Ref.  338). 


Theoretical  calculations  of  stable  pool  film  boiling  were  carried  out 
using  the  following  form  which  is  applicable  to  the  convective  contribu¬ 
tion  for  several  geometries: 


(q/A)  L 

eq 


hy”  Py  (Py  "  Pj^)  g  ^  (^  T ) 


11/4 
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The  results  are  plotted  for  each  of  three  pressures  for  each  of  the  three 

cryogenic  fluids  in  Fig,  18,  19,  and  20.  Values  of  X  were  computed 

from  Eq.  43.  The  curves  are  applicable  to  Eq .  40,  48,  and  50  i^rhere  the 

values  of  C_  and  L  are: 

5  eq 


Ge ometry 

Equation 

^5 

L 

eq 

Horizontal  Cylinder 

40 

0.62 

D 

Horizontal  Plate 
(facing  up) 

48 

0.425 

B 

(Eq.  49) 

Vertical  Surface 

50 

0.80 

L 

FORCED  CONVECTION  SUBCOOLED  NUCLEATE  BOILING 

The  previous  discussions  of  boiling  theory  were  restricted  to  the  case 
of  saturated  pool  boiling.  Boiling  can  also  occur  in  a  subcooled  liquid 
if  it  is  brought  into  contact  with  a  surface  with  a  temperature  suffici¬ 
ently  above  che  saturation  temperature.  This  is  sometimes  referred  to 
as  surface  or  local  boiling  because  the  vapor  bubbles  formed  are  subse¬ 
quently  condensed  by  the  subcooled  liquid  and  no  net  generation  of  vapor 
occurs.  The  case  where  the  subcooled  liquid  is  forced  to  flow  past  the 
heated  surface  is  of  great  practical  interest.  This  is  referred  to  as 
forced  convection  subcooled  (or  surface)  boiling.  As  with  pool  boiling, 
a  nucleate  regime  exists  which  is  limited  to  a  maximum  heat  flux  and 
transition  to  film  boiling  can  occur.  The  use  of  forced  convection  per¬ 
mits  the  attainment  of  much  higher  heat  fluxes  than  can  be  achieved  under 
poo]  conditions.  The  nucleate  regime  is  discussed  in  this  section. 
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Heat  Flux,  ,  Btu/hr-ft 


FIG.  18 


STABLE  POOL  FILM  BOILING  OF  HYDROGEN 


PIG.  19 


STABLE  POOL  FILM  BOILING  OP  NITROGEN 


Dashed  Lines  Represent  Data  of  Banchero,  ©t*  al,.  Ref 


Curves  of  heat  flux  vs  A  T  are  also  used  to  correlate  forced  convection 

suhcooled  nucleate  boiling,  McAdams  (Ref.  202)  shows  that  the  proper 

temperature  difference  to  use  is  the  wall  superheat,  At  =  T  -  T  rather 

w  s 

than  the  total  over-all  temperature  difference  between  the  wall  and  the 
subcooled  liquid  (T^  -  T^).  A  typical  curve  for  forced  convection  sub¬ 
cooled  nucleate  boiling  is  shown  in  Fig.  21,  As  discussed  by  Bergles 
and  Rohsenow  (Ref.  339),  three  distinct  regions  exist.  At  low  wall 
superheat-,  boiling  is  essentially  absent  and  the  heat  transfer  is 
governed  by  forced  convection  (Curve  AB).  At  moderate  wall  superheat, 
the  heat  transfer  is  determined  by  the  combined  effects  of  forced  con¬ 
vection  and  surface  boiling  (Curve  BC),  At  higher  wall  superheat,  the 
effect  of  forced  convection  seems  to  disappear  and  the  heat  transfer  is 
governed  by  fully  developed  nucleate  boiling  (Curve  CD).  At  some  heat 
flux  above  the  point  D,  the  maximum  heat  flux  occurs.  Several  methods 
for  predicting  the  curve  of  Fig.  21  have  been  developed. 

Rohsenow  (Ref,  30l)  recommended  an  additive  heat  flux  method  based  on 
ordinary  forced  convection  without  phase  change  and  pool  nucleate  boiling 
as  given  by  the  equation: 


(q/A)  =  (q/A)pj,  f  (q/A)g 

where  (qA)^  is  assumed  equal  to  (qA)T,T> 

D  '  PB 

(q/A)pc  .  (T„  -  T^) 


(57) 


(58) 


For  a  typical  case  of  turbulent  flow  inside  straight,  smooth,  circular 
tubes  at  moderate  temperature  differences,  the  McAdams  equation  (Ref.  202) 
gives  hp^  as: 


(59) 
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The  pool  nucleate  hoiling  contribution,  (q/A)pg,  is  calculated  from 
Rohsenow's  theory  ,  Eq.  2.  In  general,  as  shown  in  Fig.  21,  (q/A)pp  = 
(q/A)B  is  negligible  at  low  A  T,  but  becomes  dominant  at  higher  A  T. 

Nishikava  and  Yamagata  (Ref.  310)  also  recommended  Eq.  57;(q/A)pg  given 
by  a  modified  form  of  their  equation  as  presented  in  Table  6.  For  sub¬ 
cooled  forced  convection  nucleate  boiling,  they  found  it  necessary  to 
account  for  a  reduced  bubble  agitation  effect.  Accordingly,  in 
Table  6  was  modified  to: 


C 


2 


512 

10 


-  Py)  ^s  g 

B  X 


(60) 


Forster  and  Greif  (Ref.  3^5) disagreed  with  the  additive  method  given  by 

Eq.  57.  They  recommended  simply  using  their  Eq.  5  in  Table  5  for  pool 

nucleate  boiling  as  long  as  (q/A)pp  is  at  least  about  1.4  times  (q/A)p, 

for  a  given  (T  -  T  ).  They  gave  no  method  for  obtaining  the  BC  portion 
w^  s 

of  the  curve  in  Fig.  21.  As  stated  by  Zuber  and  Fried  (Ref.  302), 
Kutateladze  (Ref.  303)  and  Michenko  (Ref.  307)  also  recommended  just  cal¬ 
culating  (q/A)  on  the  basis  of  pool  nucleate  boiling  equations  alone. 


More  recently,  Kutateladze  (Ref.  3^l)  presented  the  following  expression 
for  the  entire  boiling  curve: 


q/A  =  (q/A) 


FC 

From  this  equation  it  is  seen  that: 
When  (q/A)pg/ 

When  (q/A)pB, 

When  (q/A)pg 


(q/A)pB 

2 

"  1  + 

1/2 


/(q/A)p^  =0.1, 

q/A  =1.005  (q/A)pp 

/(q/A)pc  =  10. 

q/A  =1.005  (q/A)pB 

;  = 

q/A  =  1.414  (q/A)pB 

(61) 
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This  theory  permits  the  prediction  of  the  entire  boiling  curve  ABCD  in 
Figo  21.  One  possible  method  for  turbulent  flow  in  straight,  smooth, 
circular  tubes  would  be  to  use  Eq.  6l  in  conjunction  with  Eq.  4  and  59o 
Equations  for  (q/A)^^  with  other  configurations  are  discussed  by  McAdams 
(Ref.  202). 


Levy  (Ref.  306)  extended  his  method  for  pool  nucleate  boiling  given  in 
Table  6  to  subcooled  forced  convection  nucleate  boiling  by  including  the 
sensible  heat  effect  giving: 


(Py  >■) 


(62) 


Presumably,  the  resulting  expression  for  (q/A)pg  would  have  to  be  coupled 
with  a  suitable  equation  for  (q/A)^^. 


The  most  recent  attempt  at  developing  a  theory  for  the  entire  boiling 
curve  in  Fig.  21  was  by  Bergles  and  Rohsenow  (Ref.  339).  They  suggested 
a  modification  of  Kutateladze ' s  expression  (Eq.  6l)  to  a  more  proper 
incipient  boiling  region  in  the  vicinity  of  point  B  in  Fig,  21.  Their 
recommended  equation  was: 

q/A  = 

The  term  (q/A)^^  as  shown  in  Fig.  21  is  the  boiling  contribution  to  the 
total  heat  at  incipient  boiling.  Bergles  and  Rohsenow  showed  that  the 
(q/A)g  contribution  as  obtained  under  flow’  conditions  may  not  be  the 
same  as  (q/A)pB  for  just  pool  boiling  on  the  identical  surface. 

In  summary,  the  preferred  type  expression  for  forced  convection  subcooled 
nucleate  boiling  includes  contributions  due  both  to  ordinary  forced  con¬ 
vection  without  boiling  and  to  the  boiling  effect.  As  is  discussed  in 
the  next  main  section,  little  experimental  data  on  cryogenics  exists  for 
this  type  of  boiling. 


1  + 


(q/A) 


B 

FC 


1/2 


(63) 
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FORCED  CONVECTION  SUBCOOLED  NUCLEATE  BOILING 
MAXIMUM  HEAT  FLUX 


Forced  convection  nucleate  boiling  is  also  limited  to  a  maximum  heat 
flux  and  corresponding  maximum  A  T.  For  subcooled  conditions,  a  number 
of  theories  for  estimating  this  maximum  have  been  proposed.  Some  of 
these  theories  have  been  applied  to  noncryogenic  bulk  boiling  under  lim¬ 
ited  conditions. 


Griffith  (Ref.  31?)  extended  his  pool  boiling  theory  (Eq.  2?)  to  sub¬ 
cooled  liquids  with  forced  convection  inside  tubes  by  multiplying  the 
factor  in  Table  7  by  the  following  empirical  expression.  This  expres¬ 
sion  corrects  for  both  subcooling  and  forced  convection. 


He  correlated  noncryogenic  data  with  Eq.  64  for  qualities  from  0  to  70 
percent. 


Zuber  and  Tribus  (Ref.  318)  extended  their  pool  boiling  theory  to  sub¬ 
cooled  liquids  and,  in  a  preliminary  fashion,  to  forced  convection  flow 
conditions  at  low  quality  by  an  additive  procedure  to  give  the  following 
equation. 
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(q/A)  =  |-  p  X  +  c  (T  -T  ) 

'max  2A  V  L  s  b  r, 

P," 


27r  f  ,  ^  .  /2Tr  \ 


Oi  ,  L 

L 


■'l 


k,  (T  -i,  ) 

X  Oi  ,  L  ^  s 
0  L 


where : 


g  (PL-Py) 


2  TT  g  a 


Og  g  (p  -p^. 


Pl  ’  A- 


O'  g 

X  =  2?T  .  ■‘?-— N 

0  s 


Kutateladze  (Ref.  342)  proposed  the  folloving  expression  for  flow  in 
tubes  to  include  the  effects  of  subcooling  and  forced  convection,  pre 
sumably  for  low-quality  conditions: 


(q/A), 


- 7 - X  =  0.083  V  - - 

"s  Sc  -  Af)  V  ®c 


1  .  0.057  ^ 


'L  (’/  - 

- X - 


where : 


g  (Pl  -  py 
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Gambill  (Ref.  345)  proposed  an  additive  method  similar  to  Griffith,  and 
Zuher  and  Trihus  where  he  combined  a  Kutateladze-type  expression  for 
the  boiling  contribution  and  a  McAdams-type  expression  for  the  forced 
convection  contribution  to  obtain  the  following  equation  for  low-quality 
conditions : 


Typically: 


X 


S 

24 


X'  =  0.023 

m  =  0.8 

n  =  0.4 


Gambill  conducted  an  extensive  comparison  of  Eq.  69  with  available 
experimental  data  for  the  nucleate  boiling  maximum  heat  flux  with  flow¬ 
ing,  subcooled,  wetting  liquids  in  the  absence  of  net  vaporization.  No 
cryogenic  data  were  included.  For  flow  through  tubes,  he  obtained  an 
average  deviation  of  about  19  percent  with  a  maximum  deviation  of  about 
80  percent. 


As  with  ordinary  fluids,  it  would  be  expected  that  the  forced  convection 
subcooled  nucleate  boiling  maximum  heat  flux  for  a  given  configuration- 
fluid  combination  would  depend  on  the  pressure,  degree  of  subcooling, 
and  flow  velocity  as  typically  given  by  Gambill *s  equation.  The  flow 
velocity  effect,  in  particular,  is  important  in  permitting  the  absorption 
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of  high  heat  fluxes.  A  graph  for  estimating  the  contribution  of  forced 
convection  to  nucleate  boiling,  including  the  maximum  condition,  is 
given  in  Fig.  22,  where  the  generalized  heat  transfer  coefficient  based 
on  Eq.  59  is  plotted  as  a  function  of  liquid  bulk  temperature  for  all 
three  cryogenic  fluids.  That  is: 


hpc  D 


0,2 


,0.8 


0.023 


,  0,6  0,8  p  0,4 

■•l  '’l  \ 


0.4 


=  <Pi\) 


(70) 


where  all  properties  are  evaluated  at  the  bulk  temperature. 

The  meager  cryogenic  data  available  for  the  maximum  heat  flux  are  dis¬ 
cussed  in  the  next  main  section. 
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PIG.  22  GENERALIZE®  FORCE®  CONVECTION  NON-BOILING  HEAT  TRANSFER  COEFFICIENT 
AT  SATURATE®  CONDITIONS 


FORCED  CONVECTION  SATURATED  BOILING 


In  the  two  previous  sections,  the  discussion  was  limited  to  conditions 
of  no  net  vaporization.  If  the  fluid  is  initially  saturated  or  if  suf¬ 
ficient  heat  transfer  to  the  bulk  fluid  takes  place  to  bring  it  to  the 
saturation  point,  then  a  net  change  of  phase  takes  place.  This  may  be 
referred  to  as  saturated  or  bulk  boiling. 

Some  effort  has  been  devoted  to  the  problem  of  bulk  boiling.  Unfortunately 
hovever,  no  general  unders'tanding  has  yet  been  obtained.  In  addition  to 
dependence  upon  numerous  variables  through  different  possible  modes  of 
boiling  (nucleate,  transition,  and  film)  the  pattern  of  the  two— phase 
flow  profoundly  affects  the  heat  transfer  rate.  Chen  (Ref.  3^6)  has 
illustrated  the  flow  patterns  possible  as  heat  transfer  with  net  vapori¬ 
zation  takes  place.  These  flow  patterns,  as  shown  in  Fig*  23,  include 
conditions  of  nucleate  boiling,  slug  flow,  annular  and  mist  flow,  and 
liquid  deficiency  at  the  wall. 

The  qTialitative  effects  of  fluid  quality  upon  the  heat  transfer  coef¬ 
ficient  have  been  discussed  in  considerable  detail  by  Zuber  and  Fried 
(Ref.  302)  and  Collier  (Ref.  34?)  vho  divide  the  bulk  boiling  phenomena 
into  three  distinct  regions  of  boiling.  Region  1  is  the  region  of 
nucleate  boiling  in  which  the  vapor  is  formed  at  the  heating  surface  in 
the  form  of  bubbles  that  grow,  detach,  and  become  dispersed  in  the  liquid,  or 
recondensed,  depending  upon  the  bulk  conditions  of  the  fluid.  It  is 
believed  that  in  this  region  the  heat  transfer  is  governed  by  the  mecha¬ 
nisms  of  nucleate  boiling.  Consequently,  the  equations  for  heat  transfer 
in  this  region  are  similar  to  those  for  pool  nucleate  boiling  with  an 
additional  contribution  included  to  account  for  bulk  velocity  and  its 
effect  upon  bubble  motion  as  discussed  above  and  given  by  equations  such 
as  Eq*  61  and  62. 
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Region  2  is  the  forced  convection  region  in  vhich  the  sing  and  annular 
and  mist  regimes  occur  with  net  vaporization*  Here  the  heat  transfer 
coefficient  increases  with  quality.  As  the  quality  increases  along  the 
duct,  the  increased  velocity  of  the  two-phase  mixture  may  suppress  the 
nucleate  boiling  process*  Beyond  this  point,  heat  transfer  becomes 
governed  solely  by  the  forced  convection  process.  The  actual  flow  pattern 
of  the  two-phase  mixture  is  not  exactly  knovTi*  One  school  of  thought  is 
that  only  a  liquid  film  is  in  contact  with  the  surface  (Ref,  3^iB)e 
Another  is  that  the  liquid  is  continuously  being  deposited  upon  and  re¬ 
entrained  from  the  heating  surface  (Ref.  yi9) * 

Region  3  is  the  liquid  deficient  region,  in  which  the  heat  transfer  coef¬ 
ficient  decreases  with  increasing  quality*  In  this  region,  a  liquid 
film  no  longer  wets  the  surface*  Heat  transfer  is  no  longer  due  to  a 
highly  conducting  liquid  film,  but  rather  to  a  poorly  conducting  gas. 

As  a  consequence,  the  heat  transfer  coefficient  decreases  sharply* 

An  excellent  discussion  of  recent  attempts  to  correlate  cryogenic  bulk 
boiling  is  given  by  Zuber  and  Fried  (Ref*  302)  and  is  not  repeated  here. 

In  a  very  recent  paper,  Chen  (Ref*  346)  made  a  major  contribution  to 
bulk  boiling  theory*  He  distinguished  between  the  various  flow  regimes 
for  convective  flow  of  saturated  fluids:  (l)  subcooled  nucleate 
boiling,  (2)  slug  flow,  (3)  annular  or  annular-mist  flow,  and  (4)  liquid 
deficient  conditions  as  discussed  above  and  shown  in  Fig.  23.  For 
vertical,  axial,  stable  flow  at  heat  fluxes  below  the  critical,  he 
developed  the  following  additive  relation  for  the  (no  slug)  annular  or 
annular-mist  flow  regime  which  extends  from  about  1-  to  70-percent  quality. 
The  relation  is  based  upon  the  McAdams  forced  convection  relation 
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(Eq.  59)  corrected  for  the  presence  of  vapor  by  the  Martinelli  two-phase 
parameter  (X^^)  and  the  Forster— Zuber  pool  boiling  expression  (Eq.  4) 
corrected  for  suppression  of  bubble  growth  due  to  flow.  Thus: 


q/A  -  ^  ^-Z  ^ 


where 


hpp  is  given  by  Eq.  59 

hp,_2  is  given  by  the  following  rearrangement  of  Eq.  4 


0.79  p  0.45  ..  0.49  0.25 


^-Z  " 


=  0.00122 


(A  (a 


0.5  „  0.29  ^0.24  ^  0.24 


F  is  a  two-phase  correction  as  a  function  of  X, ,  given  in  Fig.  24.  It 

XX 

accounts  for  an  increase  in  convective  turbulence  due  to  the  presence  of 


vapor 


^tt  =  Martinelli  two-phase  parameter 


S  is  a  suppression  factor  as  a  function  of  Re^^F^*  ^  given  in  Fig.  25. 
It  accoimts  for  the  suppression  of  bubble  growth  due  to  flow. 


Eep  = 


D  G  1  -  X 
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FIG.  25  SUPPRESSION  FACTOR,  S,  OP  CHEN 


Chen  tested  Eq.  71  against  existing  data  for  vater,  methanol,  cyclohexane, 
pentane,  benzene,  and  heptane  and  obtained  a  remarkable  correlation  with 
an  average  deviation  of  only  ±11  percent  for  a  quality  range  of  1  to  71 
percent.  More  theoretical  development  is  required  to  extend  Chen's 
theory  to  other  flow  regimes  and  geometric  configurations. 


Also  of  great  interest  in  forced  convection  saturated  boiling  is  the 
maximum  heat  flux.  No  generally  applicable  equations  have  yet  been 
developed.  The  effects  of  the  stream  quality  and  the  transition  from 
the  forced  convection  region  to  the  liquid  deficient  regime  add  addi¬ 
tional  complications.  An  initial  theoretical  study  of  this  problem 
has  recently  been  presented  by  Chang  (Ref.  322), 


FORCED  CONVECTION  FILM  BOILING 

General  theoretical  equations  for  forced  convection  film  boiling  both 
under  subcooled  and  saturated  flow  conditions  have  not  yet  been  developed. 
A  few  theories  for  specific  geometries  have  been  presented. 


Bromley  et  al.  (Ref.  350)  presented  tlie  following  equations  for  saturated 
forced  convection  film  boiling: 


For 


U 


<  1.0 


h  =  0.62 


V  (Pl  -  Pv)  ^L 

D  (AT)  Pr^, 


1/4 


4  r 


(75) 
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For 


>  2.0 


h  =  2.7 


1/2 


vhere 

X’  is  given  by  Eq.  42 


(76) 


and  values  for 

U' 

1.0  s  S  2.0 
Vg  D 

may  be  estimated  from  empirical  plots  in  Ref.  350. 

Motte  and  Bromley  (Ref.  35l)  later  presented  the  following  relation  for 
subcooled  film  boiling: 

,  r  D(Al)ll/"  7.29  "V 'V 
[u'  X'  I  h  D  (At) 

.  C  (A  Cj^  p.  //  ^  (C  '  (77) 

where 

C  =;  constant 

=  incident  velocity  of  liquid  on  tube 

U*  *  =  velocity  of  liquid  in  flow  passage  where  level  of 

turbulence  is  determined 
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Equations  75  to  77  are  specific  for  the  case  of  upward  flow  of  the  fluid 
outside  of  a  horizontal  heated  tube,  and  have  been  used  to  obtain  correla¬ 
tions  with  data  from  fluids  such  as  n-heptane,  carbon  tetrachloride,  and 
ethyl  alcohol.  These  equations  should  not  be  considered  as  exact  rela¬ 
tions,  and  application  to  cryogenics  has  yet  to  bo  proved. 

Cess  and  Sparrow  (Ref.  552)  obtained  a  numerical  solution  to  the  problem 
of  saturated  laminar  film  boiling  with  two-phase  flow  over  a  flat  plate 
by  a  boundary  layer  analysis.  A  closed— form  solution  ■which  approximates 
the  exact  numerical  solution  to  about  1  percent  is  as  follows: 


Nu 


=  0.5 


(P  M)l  Cy 

(P  M)y  ^ 


-1/2 


(78) 


No  experimental  data  were  presented  for  comparison,  and  the  assumptions 
necessary  for  the  solution  limit  the  applicability.  Cess  and  Sparrow 
(Ref.  353)  also  extended  the  above  case  to  subcooled  conditions  but  only 
obtained  numerical  solutions. 
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EXPERIMENTAL  DATA 


In  the  previous  section  the  various  theoretical  correlations  for  predict¬ 
ing  boiling  heat  transfer  were  reviewed.  In  most  cases,  specific  theo¬ 
retical  calculations  for  the  three  cryogenic  fluids  of  interest  were 
carried  out  and  graphical  results  were  presented.  In  this  section,  the 
available  experimental  cryogenic  boiling  data  from  the  74  references 
listed  in  the  Experimental  Data  Bibliography  are  discussed  and  compared 
with  the  theoretical  correlations. 


POOL  BOILING 

Experimental  data  for  pool  boiling  of  N^,  and  0^  are  presented  in 
Fig.  26,  27,  and  28,  respectively,  as  heat  flux  vs  A  T.  The  numbers  in 
parentheses  are  the  Bibliography  reference  numbers.  The  operating  pres¬ 
sures  are  also  stated.  For  a  given  reference,  additional  labeling  of 
the  curves  is  given  where  more  than  one  geometric  configuration  was 
tested.  Both  nucleate  and  stable  film  boiling  data  are  included.  In 
some  cases,  the  lower  end  of  the  nucleate  boiling  curves  correspond  to 
the  natural  convection  region.  Observed  values  of  the  pool  nucleate 
boiling  maximum  heat  flux  are  indicated  by  a  circle  symbol  at  the  upper 
terminus  of  the  corresponding  nucleate  boiling  curve.  Observed  values 
of  the  pool  film  boiling  minimum  heat  flux,  where  available,  are  indi¬ 
cated  by  a  square  symbol  at  the  lower  terminus  of  the  corresponding  film 
boiling  curve.  In  general,  the  curves  plotted  in  Fig.  26,  27,  and  28 
are  the  respective  authors’  correlating  lines  for  their  own  data.  Data 
covering  unusual  conditions  such  as  greased  or  roughened  surfaces  and 
zero-gravity  are  not  plotted  in  these  figures  but  are  discussed  in  the 
Conclusions  section  below. 
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Figure  26.  Experimental  Pool  Boiling  Data  for  Hydrogen 


omley 


Figure  27.  Frperimental  Pool  Boiling  Data  for  Nitrogen 


(46)  228 


Figure  28.  Experimental  Pool  Boiling  Data  for  Oxygen 


Hydrogen 


The  experimental  data  for  hydrogen,  as  plotted  in  Fig.  26,  cover  the 
following  references:  1,  15,  l6,  22,  23,  48,  54,  55,  56,  58,  62,  68, 

70,  and  73. 

As  shown,  the  nucleate  boiling  data  of  various  investigators  for  a 
variety  of  geometric  shapes  at  a  pressure  of  about  14.7  psia  (ll.5  to 
14.7)  cover  a  considerable  span  of  A  T.  For  example,  a  heat  flux  of 
1000  Btu/hr-sq  ft  was  obtained  at  a  A  T  of  about  0.5  R  by  Mulford,  et 
al.  (Ref.  54)  for  a  horizontal  cylinder,  and  at  a  A  T  of  about  5.0  R  by 
Sherley  (Ref.  62)  for  a  horizontal  flat  plate.  This  span  is  plotted  as 
dashed  lines  for  comparison  to  the  theoretical  predictions  in  Fig.  4, 

The  span  encompasses  all  the  theoretical  predictions  except  that  of 
Nishikawa  and  Yamagata,  which  appears  to  be  in  poor  agreement  with  the 
other  theories.  From  the  data  of  Class,  et  al.  (Ref.  I6)  at  12.0  psia, 
the  effect  of  surface  orientation  is  negligible.  However,  surface 
orientation  is  seen  to  assume  some  importance  at  higher  pressures.  The 
data  of  Class,  et  al.  (Ref.  I6),  Roubeau  (Ref.  58),  and  Weil,  et  al. 

(Ref.  73)  indicate  that  for  a  given  At,  the  effect  of  increasing  pres¬ 
sure  is  to  increase  the  heat  flux  as  long  as  the  maximum  is  not  approached. 
This  is  in  agreement  with  the  trend  of  theoretical  prediction  shown  in 
Fig.  7  where  the  Forster  and  Zuber  equation  is  plotted  for  three  differ¬ 
ent  pressures. 

The  data  for  the  maximum  heat  flux  as  shown  in  Fig.  26  cover  a  range 
from  about  10,000  to  43,000  Btu/hr-sq  ft  depending  mainly  on  pressure. 

These  data  are  compared  to  theoretical  predictions  in  Fig.  12.  The 
data  appear  to  scatter.  At  the  highest  pressures,  the  data  all  fall 
somewhat  below  all  the  theories. 
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The  experimental  results  for  stable  film  boiling  of  hydrogen  as  showi  in 
Fig.  26  include  several  geometric  shapes  including  wires ,  cylinders,  and 
plates.  The  plate  data  cover  a  range  oi  pressure.  As  predicted  theo¬ 
retically  in  Fig.  18,  the  film  boiling  heat  flux  shown  in  Fig.  26  does 
increase  somewhat  with  increasing  pressure.  The  data  of  Mulford,  et  al. 
(Ref.  54)  for  a  horizontal  cylinder  are  included  as  a  dashed  line  in 
Fig.  18  for  comparison  to  Bromley's  theory.  Under  Mulford's  conditions, 
the  radiation  effect  was  negligible.  Somewhat  higher  experimental  values 
were  obtained  than  were  predicted  theoretically.  Class,  et  al.  (Ref.  I6) 
reported  a  similar  difference  when  comparing  their  data  with  Bromley's 
theory.  In  general,  for  a  plate,  the  data  of  Class,  et  al.  (Ref.  I6) 
show  that  the  effect  of  going  from  a  horizontal  to  a  vertical  orientation 
is  to  increase  the  heat  flux  for  a  given  A  T.  No  experimental  pool  film 
boiling  minimum  heat  flux  values  are  available  for  comparison  to  theory. 


Nitrogen 

The  experimental  data  for  pool  boiling  of  nitrogen  as  ploi:  -fi  :.3  Fig.  27 
cover  the  following  references:  2,  3,  5,  6,  7,  9,  l<^  11,  12,  13, 

14,  18,  24,  25,  26,  27,  28,  29,  30,  31,  32,  36,  37,  41,  42,  45,  46,  47, 
49,  50,  51,  52,  58,  60,  6l,  68,  69,  71,  and  72. 

As  is  evident  from  the  above  list  and  from  a  comparison  of  Fig.  26,  27, 
and  28,  the  pool  boiling  of  nitrogen  has  been  studied  more  extensively 
than  either  hydrogen  or  oxygen. 

Compared  to  hydrogen,  the  nucleate  boiling  data  of  various  investigators 
for  nitrogen  at  essentially  14.7  psia  cover  a  much  narrower  span  on  the 
log  At  scale.  At  a  heat  flux  of  1000  Btu/hr~sq  ft,  the  At  for  nitrogen 
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is  greater  than  for  hydrogen.  It  varies  from  a  value  of  about  2.5  R 
from  the  data  of  Covley  (Ref.  18)  for  a  horizontal  cylinder,  to  a  value 
of  about  6  R  from  the  data  of  Ruzicka  (Ref.  6o)  for  a  vertical  cylinder. 
This  span  is  plotted  as  dashed  lines  for  comparison  to  the  theoretical 
predictions  in  Fig.  5.  The  span  encompasses  all  the  theoretical  correla¬ 
tions  vith  the  exception  of  the  equations  of  Kutateladze,  Chang  and  Snyder, 
and  Nishikawa  and  Yamagata.  Based  on  the  data  of  Lyon  (Ref.  45),  the 
heat  flux  for  a  given  T  is  somewhat  higher  for  a  horizontal  surface 
compared  to  a  vertical  surface.  The  data  of  Lyon,  et  al.  (Ref.  46), 
Roubeau  (Ref.  58),  and  Weil,  et  al.  (Ref.  72)  indicate  that,  in  general, 
the  effect  of  increasing  pressure  is  to  increase  the  heat  flux  for  a 
given  At.  As  with  the  other  cryogenics,  this  is  in  agreement  with  the 
trend  of  the  theoretical  prediction  as  shown  in  Fig.  8  for  nitrogen 
where  the  Forster  and  Zuher  equation  is  plotted  for  three  different 
pressures . 

The  experimental  values  for  the  maximum  heat  flux  over  a  wide  range  of 
pressure  are  compared  to  theoretical  predictions  in  Fig.  13 .  The  data 
at  a  pressure  of  about  14.7  psia  vary  over  a  very  wide  range  from  about 
18,500  Btu/hr-sq  ft  from  the  data  of  Flynn,  et  al.  (Ref,  26)  for  a 
horizontal  cylinder,  to  61,000  Btu/hr-sq  ft  from  the  data  of  Frederking 
(Ref.  28)  for  a  horizontal  wire.  In  addition,  over  the  entire  range 
of  pressure,  the  data  of  Roubeau  (Ref.  58)  for  a  horizontal  plate  differ 
markedly  from  the  data  of  Lyon,  et  al.  (Ref.  46)  for  the  same  geometrical 
condition. 

Four  data  points  for  the  minimum  heat  flux  at  atmospheric  pressure  are 
compared  to  theoretical  predictions  in  Fig,  l6.  Again,  a  considerable 
spread  in  the  data  is  evident.  For  the  most  part,  however,  the  data 
are  in  the  range  of  the  theories. 
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The  experimental  data  for  stable  film  boiling  of  nitrogen  as  shown  in 
Fig.  27,  dramatically  illustrate  the  large  effect  of  cylinder  diameter 
in  going  from  the  0.350-inch"'diameter  cylinder  of  Bromley  (Ref.  6)  to 
the  13 . 4-‘micr on-diameter  wire  of  Frederking  (Ref,  27)*  The  atmospheric 
pressure  data  of  Bromley  are  compared  with  theory  in  Fig.  19.  Excellent 
agreement  is  obtained  for  AT  ^  1000  R.  At  higher  AT  values,  the  data 
deviate  because  of  the  rapidly  increasing  contribution  of  radiation  to 
the  over-all  heat  transfer. 


Oxygen 

The  pool  boiling  experimental  data  for  oxygen  as  plotted  in  Fig.  28  cover 
the  following  references:  2,  3,  4,  34,  37,  45,  46,  47,  48,  52,  and  59. 

Compared  to  the  hydrogen  data  shown  in  Fig.  26,  the  nucleate  boiling 
data  of  various  investigators  for  oxygen  at  a  pressure  of  about  14.7  psia 
cover  a  considerably  narrower  span  on  the  log  At  scale.  At  a  heat  flux 
of  1000  Btu/hr-sq  ft,  the  A  T  for  oxygen  is  greater  than  for  hydrogen  and 
just  slightly  greater  than  for  nitrogen,  as  predicted  by  all  theories. 

It  varies  from  a  value  of  about  3  R  from  the  data  of  Lyon  (Ref.  45)  for 
a  horizontal  cylinder,  to  a  value  of  about  7  R  from  the  data  of  Bochirol, 
et  al.  (Ref.  4)  for  a  horizontal  wire.  This  span  is  plotted  as  dashed 
lines  for  comparison  to  the  theoretical  predictions  in  Fig.  6.  As  with 
nitrogen,  the  span  encompasses  all  the  theoretical  correlations  with 
the  exception  of  the  equations  of  Kutateladze,  Chang  and  Snyder,  and 
Nishikawa  and  Yamagata.  Based  on  the  data  of  Lyon  (Ref.  45)  and 
Haselten,  et  al.  (Ref.  37),  the  heat  flux  for  a  given  A  T  is  somewhat 
higher  for  a  horizontal  surface  compared  to  a  vertical  surface.  The 
data  of  Lyon,  et  al.  (Ref.  46)  indicate  that,  in  general,  the  effect  of 


increasing  pressure  is  to  increase  the  heat  flux  for  a  given  A  T.  As 
with  the  other  cryogenics,  this  is  in  agreement  with  the  trend  of  the 
theoretical  prediction  as  shown  in  Fig.  9  for  oxygen  where  the  Forster 
and  Zuher  equation  is  plotted  for  three  different  pressures. 

The  experimental  data  of  Lyon,  et  al.  (Ref.  45  and  46)  for  the  maximum 
heat  flux  over  a  wide  range  of  pressure,  are  compared  to  theoretical 
predictions  in  Fig.  14.  The  data  cover  a  range  of  about  45,000  to 
150,000  Btu/hr-sq  ft.  These  values  are  considerably  higher  than  for 
hydrogen  and  somewhat  higher  than  for  nitrogen.  As  with  hydrogen  and 
nitrogen  at  high  pressures,  a  sharper  decrease  in  maximum  heat  flux  is 
experienced  than  is  predicted. 

The  experimental  results  for  stable  film  boiling  of  oxygen  as  shown  in 
Fig.  28  include  data  on  wires  and  cylinders.  As  previously  discussed, 
the  film  boiling  heat  flux  increases  somewhat  with  increasing  pressure 
as  shown  by  the  horizontal  cylinder  data  of  Banchero,  et  al.  (Ref.  5). 

The  data  for  the  two  extreme  pressures  are  shown  as  dashed  lines  in 
Fig.  20  for  comparison  to  Bromley’s  theory.  Again,  the  radiation  con¬ 
tribution  is  negligible  over  the  range  of  the  experimental  data.  The 
agreement  is  quite  good.  As  discussed  above,  however,  a  diameter  correc¬ 
tion  is  required  to  correlate  wire  data.  No  experimental  data  for  the 
minimum  heat  flux  are  available. 

SUBCOOLED  FORCED  CONVECTION  BOILING 

The  available  experimental  data  for  subcooled  forced  convection  boiling 

are  presented  in  Fig.  29.  The  heat  flux  is  plotted  against  (T  -  T,  ) 

o 

rather  than  A  T  =  T  -  T  ,  because  in  general,  T,  was  more  well  defined 

w  S  0 

in  the  experiments  than  T  .  As  discussed  above,  however,  the  use  of 

S 

A  T  is  preferred  for  generalized  correlation. 
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Figure  29*  Experimental  Subcooled  Forced  Convection  Boiling  Data 


Data  for  both  hydrogen  and  nitrogen  are  included.  No  oxygen  data  are 
available.  The  data  cover  the  following  references:  17,  20,  21,  39, 
and  44.  References  17,  39,  and  44  also  included  bulk  boiling  data  which 
are  discussed  in  the  next  section.  The  numbers  in  parentheses  refer  to 
the  Bibliography  reference  numbers.  Other  pertinent  data  such  as  vel¬ 
ocity  range,  pressure  range,  geometry,  and  orientation  are  also  included 
in  Figo  29o  Although  the  experimental  information  is  quite  meager,  the 
nonboiling,  nucleate  boiling,  transition,  and  film  boiling  regions  are 
all  included  to  at  least  some  extent „  In  addition,  some  observed  values 
of  the  forced  convection  subcooled  nucleate  boiling  maximum  heat  flux 
are  indicated  by  circles  at  or  near  the  upper  terminus  of  the  correspond¬ 
ing  nucleate  boiling  curve » 

In  most  cases,  the  curves  plotted  in  Fig,  29  are  the  authors’  correlat¬ 
ing  lines  for  their  own  data.  With  the  exception  of  Ref.  20  and  21, 
saturated  forced  convection  data  were  also  obtained  in  the  respective 
studies,  but  the  curves  shown  in  Fig.  29  represent  subcooled  data  only. 
The  saturated  or  bulk  boiling  data  are  discussed  in  the  next  section. 

The  data  of  Lewis,  et  al.  (Ref.  44)  as  shown  in  Fig.  29  are  for  the 
nucleate  boiling  of  nitrogen  flowing  upward  inside  a  vertical  circular 
tube  having  an  inside  diameter  of  0.555  inch.  The  degree  of  subcooling 
was  less  than  5  R  and  both  the  pressures  and  flow  velocities  were  rela¬ 
tively  low.  At  these  conditions,  the  contribution  of  forced  convection 
to  the  total  heat  flux  as  given  in  Eq.  58,  59,  and  6l,  and  as  estimated 
from  Fig.  22  (assuming  turbulent  flow),  is  less  than  1000  Btu/hr-sq  ft. 
This  is  much  less  than  the  lowest  heat  flux  shown  for  the  Lewis  datao 
Thus,  the  forced  convection  process  is  strongly  dominated  by  nucleate 
boiling.  As  would  be  expected,  therefore,  the  Lewis  heat  flux  curve 
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for  Ng  shown  in  Figo  29  agrees  closely  with  the  low  pressure  pool  boil¬ 
ing  curves  for  nitrogen  shown  in  Fig,  27.  Several  values  of  Lewis’ are 
also  shown  in  Fig.  29  for  the  maximum  nucleate  boiling  heat  flux.  Again, 
for  the  low  velocities  studied,  the  forced  convection  contribution  to 
the  maximum  heat  flux  as  given  by  Eq.  69  is  relatively  small.  Compari¬ 
son  of  these  maximum  heat  flux  points  with  similai*  pool  boiling  data  in 
Fig.  27  shows  fairly  good  agreement;  the  Lewis  values  appearing  to  be 
generally  somewhat  low.  Lewis,  et  al,  also  obtained  some  hydrogen  data, 
but  did  not  consider  the  temperature  measurements  to  be  of  suffici^ent 
accuracy  to  warrant  any  correlation. 


The  data  of  Dean,  et  al,  (Ref.  20  and  2l)  as  shown  in  Fig,  29  are  for 

nitrogen  flowing  through  a  vertical  annulus  with  an  equivalent  diameter 

(D  -  D.)  of  0,138  inch.  The  data  cover  the  entire  region  from  nonboil- 
0  1 

ing  to  film  boiling.  The  degree  of  subcooling  was  of  the  order  of  50 
to  60  F  at  a  pressure  of  300  psia  and  a  relatively  high  flow  velocity 
of  40  ft/sec.  For  the  nucleate  boiling  region  at  these  conditions, 
the  contribution  of  forced  convection  to  the  heat  flux  may  be  estimated 
from  Fig,  22  by  using  the  equivalent  diameter  and  modifying  the  ordinate 
as  recommended  for  an  annulus  by  Wiegand  as  given  in  McAdams  (Ref.  202). 
Thus,  the  ordinate  in  Fig.  22  becomes: 

yO.S  ^ 


where  in  this  case,  D.  =  0.01568  ft  and  D  =  0,02717  ft.  A  value  of 

1  0 

(q/A)j,p  of  the  order  of  350,000  is  accordingly  estimated  for  T^  -  T^  = 
60  R.  This  is  of  the  order  of  the  experimental  nonboiling  and  nucleate 
boiling  results  of  Dean  as  shown  in  Fig,  29*  At  this  relatively  high 
velocity,  the  contribution  of  forced  convection  to  the  total  heat  flux 
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is  very  significant.  The  nucleate  boiling  contribution  to  the  total 

heat  flux  is  difficult  to  estimate  because  of  experimental  uncertainty 

as  to  the  T  -  T  values.  However,  for  a  pressure  of  300  psia,  the  max- 
w  s 

imum  pool  nucleate  boiling  heat  flux  for  nitrogen  as  estimated  from 
Figo  13  is  of  the  order  of  100,000  Btu/hr-sq  ft.  The  contribution  due 
to  subcooling  is  estimated  from  Eq,  69  to  be  about  15,000  Btu/hr-sq  ft. 
Thus,  the  total  estimated  heat  flux  in  the  region  of  the  maximum  is  about 
465,000  Btu/hr-sq  ft.  This  is  in  fairly  good  agreement  with  the  data „ 

The  results  of  Dean  also  extend  into  the  film  boiling  region  where  no 
suitable  theoretical  correlation  is  available  for  comparison.  Compared 
to  the  pool  film  boiling  data  of  Fig.  27,  however,  it  is  evident  that  a 
significant  effect  of  forced  convection  exists.  In  summary,  comparison 
of  the  forced  convection  results  of  Dean,  et  al.  in  Fig.  29  with  pool 
boiling  results  in  Fig.  27  shows  that,  as  predicted  by  theory,  the  appli¬ 
cation  of  forced  convection  in  conjunction  with  the  boiling  process  per¬ 
mits  the  attainment  of  much  higher  heat  fluxes  than  are  possible  with 
pool  boiling  only. 

The  data  of  Core,  et  al.  (Ref.  17 )  as  shown  in  Fig.  29  are  for  hydrogen 

flowing  through  a  vertical  circular  tube  with  an  inside  diameter  of 

0,1675  incho  The  data  cover  a  region  from  nucleate  to  film  boiling.  The 

degree  of  subcooling  was  as  high  as  approximately  15  F.  The  pressure 

was  varied  from  about  JO  to  I60  psia  and  the  velocity  from  about  25  to 

60  ft/sec.  At  Coreas  operating  conditions,  the  contribution  of  forced 

convection  to  the  heat  flpx  is  estimated  from  Fig.  22  to  vary  from  about 

50,000  to  200,000  Btu/hr-sq  ft  over  the  range  of  pressure,  velocity,  and 

T  “  T,  studied.  In  general,  as  seen  in  Fig,  29,  this  range  of  heat  flux 
w  0 

values  appears  to  be  higher  than  measured  by  Core,  et  al.  for  the  nucleate 
boiling  region.  The  apparent  failure  of  Eq.  70  and  Fig.  22  in  this  case 
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to  satisfactorily  predict  (q/A)^^^  may  be  (hic  to  the  veil-known  difficulty 
of  making  estimations  for  liquids  in  the  vicinity  of  the  critical  tem¬ 
perature,  especially  when  the  wall  temperature  is  above  the  critical, 
as  was  the  case  here.  In  this  region,  physical  properties  can  change 
drastically  and  the  choice  of  the  bulk  temperature  to  evaluate  physical 
properties  is  very  uncertaiuo  Again,  it  is  difficult  to  estimate  the 
nucleate  boiling  contribution  to  the  heat  flux,  but  the  maximum  values 
for  pool  boiling  as  seen  in  Fig,  12  and  26  vary  from  about  15,000  to 
40,000  Btu/hr-sq  ft,  depending  on  the  pressure.  These  values  are  in  the 
lover  range  of  Core’s  datao  Thus,  it  appears  probable  that  both  nucleate 
boiling  and  forced  convection  were  contributing  to  the  over-all  heat  flux,, 

The  data  of  Hendricks,  et  al.  (Ref.  39)  as  sho^vai  in  Fig.  29  are  for 
hydrogen  flowing  through  a  vertical  circular  tube  with  an  inside  diameter 
of  0.313  incho  The  subcooled  data  cover  only  the  film  boiling  region. 

The  degree  of  subcooling  was  less  than  about  5  F„  The  pressure  was  var¬ 
ied  from  about  50  to  70  psia  and  the  velocity  from  about  70  to  85  ft/sec . 
Compared  to  the  pool  film  boiling  data  shoivTi  in  Fig^  26,  much  higher 
heat  fluxes  were  obtained  by  Hendricks  with  forced  convection  f or  a 
given  temperature  difference „  The  data  of  Core  in  the  film  boiling 
region  are  in  approximate  agreement  with  that  of  Hendricks. 


SATURATED  FORCED  CONVECTION  BOILING 

Experimental  data  for  the  saturated  forced  convection  boiling  of  hydrogen, 
nitrogen,  and  oxygen  are  covered  in  the  following  references:  17,  35, 

39^  43^  kk,  53,  57,  63,  66,  67,  and  74.  As  discussed  above,  this  area 
of  cryogenic  boiling  is  considerably  more  complicated  than  the  areas  of 
pool  boiling  or  subcooled  forced  convection  boiling.  Compared  to  the 
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latter,  at  least  two  additional  variables  must  be  considered  in  saturated 
forced  convection  boilings  These  are  quality  of  the  flowing  stream  and 
type  of  flow  regime.  Because  of  the  large  number  of  parameters  involved, 
and  the  general  lack  of  sufficient  study  of  each  variable  independently, 
the  data  are  not  plotted  as  heat  flux  vs  temperature  here.  Instead,  the 
ranges  of  parameters  covered  in  the  various  investigations  are  summarized 
in  Table  9  and  the  findings  are  discussed. 

Core,  et  al.  (Ref.  I7)  obtained  data  for  hydrogen  at  relatively  high 
velocities  over  a  wide  range  of  pressure  and  quality.  Their  data  covered 
a  wide  range  of  boiling  conditions.  Data  in  the  nucleate  regime  were 
obtained  at  surprisingly  high  values  of  T^  -  T^  (50  to  100  R).  Heat 
fluxes  as  high  as  3jl00,000  Btu/hr-sq  ft  were  achieved  in  the  film  boil¬ 
ing  region. 

Hendricks,  et  al.  (Ref.  35  and  39)  also  obtained  data  for  hydrogen  at 
relatively  high  velocities  over  a  wide  range  of  pressure  and  quality. 
However,  they  covered  only  the  film  boiling  region.  Heat  fluxes  as  high 
as  520,000  Btu/hr-sq  ft  were  obtained. 

Jones,  et  al,  (Ref.  43)  studied  heat  transfer  to  saturated  nitrogen  at 
relatively  low  velocities  in  a  very  long  0.43-inch-ID  horizontal  circu¬ 
lar  tube.  They  determined  the  critical  quality  corresponding  to  the 
point  at  which  the  heat  transfer  coefficient  drops  off  sharply  due,  pre¬ 
sumably  to  transition  to  the  liquid  deficient  region  shown  in  Fig,  23 
and  discussed  above.  Heat  fluxes  to  about  4500  Btu/hr-sq  ft  were  meas¬ 
ured,  The  critical  quality  varied  from  80  to  99  percent,  depending  on 
heat  flux  and  flowrate. 
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Lewis,  et  al.  (Ref.  44  and  66)  conducted  tests  with  both  hydrogen  and 
nitrogen  at  relatively  low  velocities.  They  also  concentrated  on  deter¬ 
mining  the  critical  quality  and  corresponding  maximum  heat  flux,  hut 
over  a  somewhat  wider  range  of  variables  than  Jones,  et  al.  With  hydro¬ 
gen,  they  achieved  maximum  nucleate  boiling  heat  fluxes  up  to  25,000  Btu/ 
hr-sq  ft  over  a  wide  range  of  critical  quality.  With  nitrogen,  maximum 
nucleate  boiling  heat  fluxes  to  26,000  Btu/hr-sq  ft  were  achieved  over 
a  wide  range  of  quality.  Data  in  the  film  boiling  region  downstream  of 
the  transition  point  were  obtained  at  heat  fluxes  up  to  40,000  Btu/hr- 
sq  ft  with  both  cryogenic  fluids. 

The  data  presented  by  Monroe,  et  al.  (Ref.  53)  cover  only  the  transition 
region  for  oxygen  and  nitrogen  flowing  in  a  vertical  tube  with  an  inside 
diameter  of  0.5  inch  at  very  low  velocity  and  low  pressure.  The  data 
of  Richards,  et  al.  (Ref.  57)  3nd  Sydoriak,  et  al.  (Ref.  63)  present 
data  for  the  maximum  nucleate  boiling  heat  flux  of  nitrogen  and  hydrogen 
for  flow  through  vertical  annuli  at  moderate  velocities  to  relatively 
high  velocities.  Critical  heat  fluxes  up  to  about  38,000  Btu/hr-sq  ft 
were  achieved. 

The  data  of  Wright  and  Walters  (Ref.  6?  and  74)  are  for  hydrogen  flowing 
in  a  horizontal  tube  with  an  inside  diameter  of  0.25  inch.  A  moderate 
pressure  range  was  covered  at  relatively  high  velocities.  Only 
low  exit  quality  values  were  obtained.  Under  these  limited  conditions, 
a  satisfactory  correlation  of  the  data  was  obtained  for  both  the  nucleate 
and  film  boiling  regions.  The  use  of  relatively  high  velocities  resulted 
in  heat  fluxes  as  high  as  123,000  Btu/hr-sq  ft  for  the  film  boiling  region. 
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CONCLUSIONS 


Based  on  the  experimental  information  contained  in  the  literature  sur¬ 
veyed  during  this  study,  and  the  comparisons  made  with  theoretical  pre¬ 
dictions,  as  discussed  in  the  previous  section,  the  following  conclusions 
may  he  drawn  concerning  heat  transfer  to  boiling  hydrogen,  nitrogen,  and 
oxygen o 


POOL  NUCLEATE  BOILING 

Data  at  Atmospheric  Pressure 

A  considerable  quantity  of  pool  nucleate  boiling  experimental  data  has 
been  obtained  for  all  three  cryogenic  fluids  at  atmospheric  pressure. 

In  general,  as  with  other  fluids,  the  heat  flux  increases  rapidly  with 
increasing  A  T.  The  data  cannot  be  correlated  by  a  single  line,  but 
cover  a  relatively  wide  range  of  A  T  for  a  given  heat  flux.  The  spread 
is  attributed  primarily  to  effects  of  surface  condition,  and  surface 
geometry  and  orientation.  In  general,  hydrogen  boils  at  much  lower  val¬ 
ues  of  A  T  than  nitrogen  or  oxygen.  Nitrogen  may  boil  at  slightly  lower 
values  of  A  T  than  oxygen.  This  is  consistent  with  the  trends  predicted 
by  theoryo 


Surface  Roughness  and  Orientation 

Based  on  the  experimental  data  of  Ref.  15,  45,  and  52  covering  all  three 
cryogenic  fluids,  a  profound  effect  of  surface  roughness  on  pool  nucleate 
boiling  exists.  In  general,  when  a  smooth  or  polished  surface  is  rough¬ 
ened,  the  required  At  for  a  given  nucleate  boiling  heat  flux  is  reduced. 
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Because  of  the  large  effect  of  surface  roughness,  detailed  comparisons 
and  firm  conclusions  for  different  geometrical  shapes  are  not  possible « 
The  experimental  data  of  Ref.  49  and  60  for  nitrogen  do  indicate  little 
effect  of  diameter  for  tubes  and  wires  on  the  nucleate  boiling  heat  flux, 
except  at  low  heat  flux  values  where  natural  convection  may  be  predomi¬ 
nant.  The  effect  of  surface  orientation  is  somewhat  uncertain.  The 
data  of  Ref.  15  for  indicate  a  negligible  effect  on  the  AT  for  a 
plate  in  going  from  a  vertical  to  a  horizontal  orientation.  On  the  other 
hand,  the  data  of  Ref.  37  and  45  for  nitrogen  and  oxygen  indicate  that, 
in  general,  the  required  A  T  for  a  given  nucleate  boiling  heat  flux  with 
a  cylinder  is  decreased  in  going  from  a  vertical  to  a  horizontal 
orientation. 


Comparison  of  Theories  With  Data 


In  general,  the  pool  nucleate  boiling  theories  summarized  in  Tables  5 
and  6  do  not  consider  the  nature,  geometry,  or  orientation  of  the  heater 
surface.  Comparison  of  the  theories  with  data  at  atmospheric  pressure 
indicates  that  the  correlations  of  Rohsenow,  Forster  and  Zuber,  Forster 
and  Greif,  Levy,  Michenko,  and  Laboutzov  fall  within  the  spread  of  the 
experimental  data. 
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Increasing  Pressure 


For  a  given  fluid-heater  combination,  sufficient  experimental  data  are 
available  to  show  that  the  effect  of  increasing  pressure  is  to  decrease 
the  required  A  T  for  a  given  pool  nucleate  boiling  heat  flux.  This  is 
in  agreement  with  the  trends  predicted  by  the  pool  nucleate  boiling 
theories  such  as  that  of  Forster  and  Zuber, 
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Effect  of  Gravity 


As  shown  both  by  the  experimental  data  of  Sherley  (Ref.  62)  for  hydrogen 
and  Clark,  et  al.  (Ref.  8)  for  nitrogen,  the  effect  of  both  very  low  and 
high  gravity  on  the  nucleate  boiling  heat  flux  is  practically  negligible. 
As  explained  by  Clark,  this  is  probably  due  to  the  relatively  high  bubble 
Froude  numbers  (ratio  of  inertial  to  gravitational  force)  for  cryogenic 
fluids . 


Greased  or  Coated  Surface 


The  experimental  data  of  Ref.  15  and  45  show  that  the  required  At  may 
be  increased  significantly  for  a  given  nucleate  boiling  heat  flux  if  the 
surface  is  greased  or  coated  with  a  glossy  film.  This  is  most  pronounced 
for  a  horizontal  orientation  of  the  surface. 


POOL  NUCLEATE  BOILING  MAXIMUM  HEAT  FLUX 


Effect  of  Pressure 


Some  experimental  data  on  the  pool  nucleate  boiling  maximum  heat  flux  arc 
available  for  all  three  cryogenic  fluids  over  a  wide  range  of  pressure. 

In  general,  as  with  other  fluids,  the  maximum  heat  flux  curve  goes 
through  a  peak  value  at  a  pressure  somewhat  less  than  50  percent  of  the 
critical  pressure.  For  a  given  fluid,  the  data  cannot  be  correlated  by 
a  single  line,  but  instead  cover  a  range  of  maximum  heat  flux  for  a 
given  pressure.  No  satisfactory  explanation  exists  for  the  spread  of 
the  data.  In  general,  at  a  given  reduced  pressure,  the  maximum 
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heat  flux  of  hydrogen  is  considerably  lover  than  for  nitrogen  and  oxygen. 
Somewhat  higher  maximum  heat  fluxes  can  be  achieved  with  oxygen  as  com¬ 
pared  to  nitrogen.  This  is  consistent  with  the  trends  predicted  by 
theory. 


Surface  Roughness  and  Orientation 

Based  on  the  experimental  data  of  Ref.  15,  45,  and  54  covering  all  three 
cryogenic  fluids,  the  effect  of  surface  roughness  on  the  maximum  4ie at 
flux  is  small  or  negligible.  If  the  surface  is  greased,  however,  the 
data  of  Ref.  15  for  hydrogen  show  that  the  maximum  heat  flux  can  be 
increased  by  from  50-  to  100-percent,  depending  on  surface  orientation. 
The  data  of  Ref.  45  for  nitrogen  and  oxygen  also  show  some  increase  in 
the  maximum  heat  flux  when  the  surface  is  coated  with  a  glossy  film. 

The  data  of  Ref.  45  further  show  that  a  higher  maximum  heat  flux  can  be 
achieved  with  a  vertical  surface  than  with  a  horizontal  surface.  This 
is  in  disagreement  with  the  opposite  prediction  of  the  equations  of 
Chang  (Ref.  322). 


Comparison  of  Theories  With  Data 


In  general,  the  pool  nucleate  boiling  maximum  heat  flux  theories  summar¬ 
ized  in  Table  7  do  not  consider  the  geometry  or  orientation  of  the  heater 
surface.  Comparison  of  the  theories  with  experimental  data  as  shown  in 
Fig.  12,  13,  and  14  shows  that  the  equations  of  Rohsenow  and  Griffith, 
Zuber  and  Tribus,  Kutateladze,  Borishanskii ,  Noyes,  Chang  and  Snyder,  and 
Moissis  and  Berenson  fall  within,  or  in  the  near  vicinity  of  the  data 
scatter,  except  at  high  pressures. 
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Effect  of  Gravity 


As  shown  by  the  experimental  data  of  Sherley  (Ref.  62)  for  hydrogen  and 
Clark,  et  al.  (Ref.  8)  for  nitrogen,  the  maximum  heat  flux  is  reduced 
by  about  50-percent  in  going  from  standard  gravity  to  near-zero  gravity. 


POOL  FILM  BOILING  MINIMUM  HEAT  FLUX 

Few  experimental  data  are  available  for  the  pool  film  boiling  minimum 
heat  flux.  Only  nitrogen  at  atmospheric  pressure  has  been  studied  to 
any  extent  and  the  results  appear  to  vary  widely,  depending  on  surface 
geometry  and  orientation.  The  data  do  fall  roughly  within  the  theoreti¬ 
cal  predictions.  Based  on  the  data  of  Lyon,  et  al,  (Ref.  46),  it  is 
interesting  to  note  that  as  the  critical  pressure  is  approached,  the 
transition  region  between  nucleate  and  film  boiling  seems  to  disappear 
leaving  only  an  inflection  point  in  the  q/A  vs  A  T  curve  in  place  of 
the  maximum  and  minimum  heat  flux  values.  Near  the  critical  pressure, 
stable  film  boiling  can  occur  at  very  low  values  of  A  T. 


STABLE  POOL  FILM  BOILING 

Data  at  Atmospheric  Pressure 

A  fair  quantity  of  stable  pool  film  boiling  experimental  data  have  been 
obtained  for  all  three  cryogenic  fluids  at  atmospheric  pressure.  In 
general,  as  with  other  fluids,  the  heat  flux  increases  with  increasing 
At  until  the  surface  fails.  At  high  surface  temperatures,  the  effect 
of  the  radiation  contribution  may  become  important.  In  general,  for  a  given 
A  T,  somewhat  higher  heat  fluxes  can  be  achieved  with  hydrogen  than  with 
nitrogen  or  oxygen.  This  is  consistent  with  the  trends  predicted  by 
theory. 
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Surface  Roughness  and  Orientation 


Based  on  the  experimental  data  of  Ref.  15,  and  as  predicted  theoretically, 
there  is  no  effect  of  surface  roughness  or  a  greased  surface  on  pool 
film  boiling.  The  experimental  data  of  Ref.  15  and  42  show  that  a  hori¬ 
zontal  surface  requires  a  higher  A  T  to  sustain  a  given  heat  flux  than 
a  vertical  surface.  For  horizontal  cylinders  or  wires,  the  data  of 
Ref.  3  and  27  indicate  an  increasingly  larger  effect  of  diameter  at 
values  less  than  0.1  inch  than  predicted  hy  the  theory  of  Bromley  as 
discussed  in  the  Boiling  Theory  section  above.  Frederking  (Ref.  2?),. 
in  particular,  has  achieved  heat  fluxes  of  almost  1  Btu/sec-sq  in.  with 

very  fine  wires. 


Effect  of  Pressure 

With  proper  corrections  for  small  diameter,  height  and/or  radiation,  the 
pool  film  boiling  theories  provide  a  relatively  adequate  means  of  pre¬ 
dicting  pool  film  boiling  heat  flux.  The  effect  of  pressure  as  shown 
hy  the  experimental  data  of  Ref.  3  is  to  increase  the  heat  flux;  this  is 
in  good  agreement  with  the  theory. 


Effect  of  Gravity 

The  data  of  Ref.  8  show  that  near-zero  gravity  conditions  cause  a  con¬ 
siderable  reduction  in  film  boiling  heat  flux  for  a  given  A  T,  as  pre¬ 
dicted  hy  theory. 


SUBCOOLED  FORCED  COJfVECTION  BOILING 


Compared  to  that  available  for  pool  boiling,  only  a  very  limited  amount 
of  data  are  available  for  subcooled  forced  convection  boiling.  No  known 
data  for  oxygen  even  exist.  Therefore,  few  conclusions  can  be  drawn. 

In  general,  the  results  available  do  show  that  forced  convection  can  be 
used  to  achieve  much  higher  heat  fluxes  than  are  possible  with  pool  boil¬ 
ing.  This  is  true  not  only  in  the  nucleate,  but  also  in  the  film  boiling 
region.  Theoretical  methods  such  as  Eq.  6l,  63,  and  69  show  promise  for 
allowing  approximate  predictions  of  the  nucleate  boiling  region. 


SATURATED  FORCED  CONVECTION  BOILING 

This  type  of  boiling  is  considerably  more  complex  than  the  other  types 
discussed  above.  Again,  only  a  limited  amount  of  experimental  data  are 
available.  Oxygen  has  only  been  studied  in  the  transition  region.  The 
results  show  that  at  relatively  low  heat  fluxes,  a  nucleate-type  boiling 
is  obtained  until  a  critical  heat  flux  and  a  corresponding  critical 
quality  are  reached.  Beyond  this  point,  a  film  boiling-type  region 
exists  which  may  or  may  not  correspond  to  a  liquid  deficient  condition, 
depending  on  the  quality  and  flow  velocity  of  the  flowing  stream.  As 
with  subcooled  conditions,  very  high  heat  fluxes  can  be  attained  with 
saturated  forced  convection  boiling. 
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RECOMMENDATIONS 


On  the  basis  of  the  results  of  the  literature  survey  presented  above, 
the  following  recommendations  are  made  for  additional  analytical  and 
experimental  studies  that  are  required  to  provide  an  adequate  knowledge 
of  heat  transfer  to  boiling  cryogenics. 

1.  One  of  the  most  important  variables  in  pool  nucleate  boiling  is 
the  nature  of  the  surface.  Additional  experiments  with  all 
three  cryogenic  fluids  are  required  with  surfaces  of  known  dis¬ 
tribution  and  geometry  of  bubble  nucleation  sites.  The  data 
should  be  compared  to  theories  such  as  that  of  Lienhard  (Ref.  313). 

2.  Additional  data  on  the  maximum  pool  nucleate  boiling  heat  flux 
for  all  three  fluids  over  a  wide  range  of  pressure,  geometry, 
and  orientation  are  required  to  determine  the  reasons  for  the 
apparent  disagreement  among  the  meager  data  thus  far  available. 

This  discrepancy  has  also  been  discussed  recently  by  Squire 
(Ref.  207). 

3.  Experimental  data  on  the  pool  film  boiling  minimum  heat  flux 
are  practically  nonexistent  and  are  needed  for  all  three  fluids 
over  a  wide  range  of  pressure  and  conditions  of  geometry  and 
orientation  for  adequate  comparison  with  available  theory, 

4o  Subcooled  forced  convection  boiling  studies  are  greatly  needed. 

They  should  be  conducted  with  all  three  cryogenic  fluids  over 
a  wide  range  of  conditions  of  pressure,  velocity,  and  heat  flux 
in  both  straight  and  curved  tubes.  Data  covering  the  entire 
region  from  nonboiling  through  nucleate,  transition,  and  to 
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stable  film  boiling  should  be  obtained.  Supplementary  pool 
boiling  measurements  should  be  made  with  the  same  materials  to 
provide  basic  information  to  test  the  adequacy  of  available 
theories  in  the  nucleate  boiling  region.  The  importance  of  this 
has  recently  been  shown  by  Bergles  and  Rohsenow  (Ref.  339). 

These  studies  can  be  conveniently  made  in  e lectrically  heated 
tubes  with  a  pressurized  flow  system.  The  concurrent  develop¬ 
ment  of  a  theory  for  the  subcooled  film  boiling  region  is 
necessary*  ^ 

The  need  for  more  detailed  and  comprehensive  experimental  and 
analytical  studies  of  saturated  forced  convection  boiling  with 
all  three  cryogenic  fluids  cannot  be  overemphasized.  Data 
should  be  obtained  in  the  same  manner  and  over  the  same  ranges 
of  variables  as  discussed  in  Item  4  above.  Local  axial  condi¬ 
tions  must  be  carefully  determined.  Both  partial  and  complete 
vaporization  of  the  fluids  should  be  examined.  Chen's  theory 
(Ref.  346)  should  be  tested  to  determine  its  adequacy  for  cryo¬ 
genic  fluids  up  to  the  critical  quality.  Data  in  the  liquid 
deficient  region  should  be  compared  to  ordinary  forced  convec¬ 
tion,  nonboiling  correlations. 


NOMENCLATURE 


C 

D 

f 

F 

g 

G 


h 

J 

k 

L* 

L 

o 

Nu 

P 

A  P 


Pr 

q/A 

Ra 


ratio  of  local  to  sea-level  gravitational  acceleration 

specific  heat  at  constant  pressure 

diameter 

frequency 

heating  surface  condition  factor 
Reynolds  number  factor  of  Chen 
acceleration  due  to  gravity 

Q 

constant  =  32.17  Ibm-ft/lbf-sec 
mass  velocity 

characteristic  mass  velocity  for  pool  boiling 
heat  transfer  coefficient 

mechanical  equivalent  of  heat  =  778*2  ft-lbf/Stu 

thermal  conductivity 

characteristic  length  for  boiling 

transition  height 

Nusselt  number 

pressure 

pressure  difference  corresponding  to  A  T 
Prandtl  number 
heat  flux 
Rayleigh  number 


i 
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Re  =  Reynolds  number 

S  =  suppression  factor  of  Chen 

St  =  Stanton  number 

T  =  temperature 

At  =  wall  superheat  =  T  -  T 

w  s 

u,  U  =  velocity 

V  =  velocity 

X  =  weight  fraction  vapor  (quality) 

X. ,  =  Martinelli  two-phase  parameter 

't't 

y  =  height 

z  =  weight  fraction  liquid 

Greek  Letters: 

Ot  =  thermal  diffusivity 

Oe’  =  absorptivity 

=  coefficient  of  volumetric  expansion 
A  =  difference 

C  =  emissivity 

X  =  heat  of  vaporization 

fX  =  viscosity 

IT  =  constant  = 

p  =  density 

a  =  surface  tension 

Cy*  =  Stefan-Boltzmann  constant 


118 


Subscripts : 


atm 

b 

B 

Bi 

c 

f 

FC 

L 

PB 

r 

s 

"L 

V 

W 


atmospheric 
bulk  condition 
boiling 

incipient  boiling 
convection 
film  condition 
forced  convection 
liquid 

pool  boiling 
radiation 
saturation 
saturated  liquid 
vapor 

wall  condition 
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Commission  I,  p  49-53 

N^:  A1-B1-C1-E2 

N^:  A4-B1-C1-E2  ^ 

H^:  A1-B1-C1-E2 
H^:  A4-B1-C1-E2 


*59 »  Rushton,  Ja  Ho  (data  of  Wo  Fo  Giaque) 

“Improved  Equipment  for  Oxygen  Production” 

NDRCj  Division  11,  Summary  Technical  Report,  Volume  1,  Chapter  7, 
Washington,  Do  Co  (l946) 

0^:  Al-(lncomplete  description)  X 

*60.  Ruzicka,  Jo 

”Heat  Transfer  to  Boiling  Nitrogen” 

Bull.  Inst.  Interne  froid.  Annexe  1958-1,  3^3-9 


N 


2‘ 


2' 


N. 


A1-B1-C5,  C9-E2 
A2-B1-C5-E2 
A3-B1-C5-E2 
A4-B1-C9-E2 


XX 


*61 


Ruzicka,  J. 

“Heat  Transfer  to  Boiling  Nitrogen" 
Technika  Hlubokych  Teplot 

Nakladatelstvi  Ceskoslovenske  Akademic  Ved 
Praha  1961,  pp  109-152 

N^:  A1-B1-.C5,  C9-E3 

N^:  A2-B1-C5-E3 
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SUBJECT  CROSS  REFERENCES 

(Numbers  Refer  to  Experimental  Data  Bibliograpl^) 


BOILING  REGIMES  AND  FLOW  CONDITIONS 

Hydrogen  Pool  Nucleate  Boiling: 

1,  15,  16,  22,  23,  48,  54,  55,  56,  58,  62,  64,  68,  70,  73 

Nitrogen  Pool  Nucleate  Boiling: 

3,  4,  8,  9,  10,  11,  12,  13,  14,  18,  24,  25,  26,  28,  33,  37,  45, 
46,  47,  49,  50,  51,  52,  58,  60,  6l,  65,  68,  69,  71,  72 

Oxygen  Pool  Nucleate  Boiling: 

3,  4,  34,  37,  45,  46,  47,  48,  52,  55,  59,  65 

Hydrogen  Subcooled  Forced  Convection  Nucleate  Boiling: 

17,  44 

Nitrogen  Subcooled  Forced  Convection  Nucleate  Boiling: 

20,  21,  44 

Hydrogen  Saturated  Forced  Convection  Nucleate  Boiling: 

17,  44,  57,  63,  67,  74 

Nitrogen  Saturated  Forced  Convection  Nucleate  Boiling: 

43,  44,  57,  63 

Hydrogen  Maximum  Pool  Nucleate  Boiling  Heat  Flux: 

1,  15,  54,  56,  58,  62 

Nitrogen  Maximum  Pool  Nucleate  Boiling  Heat  Flux: 

8,  9,  10,  11,  12,  13,  14,  18,  26,  28,  45,  46,  47,  50,  51,  54,  58, 
60,  61 

Oxygen  Maximum  Pool  Nucleate  Boiling  Heat  Flux: 

45,  46,  47 

Hydrogen  Maximum  Forced  Convection  Nucleate  Boiling  Heat  Flux: 

44,  66 

Nitrogen  Maximum  Forced  Convection  Nucleate  Boiling  Heat  Flux: 

20,  44,  66 
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Hydrogen  Transition  Boiling: 

17 

Nitrogen  Transition  Boiling: 

8,  9,  10,  11,  12,  13,  14,  18,  20,  26,  50,  51,  53,  60,  6l 
Oxygen  Transition  Boiling: 

53 

Hydrogen  Pool  Film  Boiling: 

15,  16,  54,  68,  70 

Nitrogen  Pool  Film  Boiling: 

3,  4,  5,  6,  7,  8,  9,  10,  11,  13,  14,  18,  26,  27,  29,  30,  31,  32, 
36,  41,  42,  46,  47,  50,  51,  60,  6l,  68,  69 

Oxygen  Pool  Film  Boiling: 

2,  3,  4,  34,  46,  47 

Hydrogen  Forced  Convection  Subcooled  Film  Boiling: 

17,  35,  39 

Nitrogen  Forced  Convection  Subcooled  Film  Boiling: 

20,  21 

Ifydrogen  Forced  Convection  Saturated  Film  Boiling: 

17,  35,  39,  40,  44,  66,  67,  74 

Nitrogen  Forced  Convection  Saturated  Film  Boiling: 

44,  66 

GEOMETRY  CONDITIONS 

Hydrogen  Boiling  above  Horizontal  Plates: 

1,  15,  16,  56,  58,  62 

Nitrogen  Boiling  above  Horizontal  Plates: 

46,  47,  58 

Oxygen  Boiling  above  Horizontal  Plates: 

46,  47 

Hydrogen  Boiling  on  Vertical  Plates: 

15,  16 

Nitrogen  Boiling  on  Vertical  Plates 
24,  25 
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Hydrogen  Boiling  Inside  Vertical  Tubes: 

17,  35,  39,  40,  44,  57,  63,  66 

Nitrogen  Boiling  Inside  Vertical  Tubes: 

20,  21,  38,  44,  53,  57,  63,  65,  66 

Oxygen  Boiling  Inside  Vertical  Tubes: 

38,  53,  65 

Hydrogen  Boiling  Outside  Vertical  Tubes: 

22,  23,  64 

Nitrogen  Boiling  Outside  Vertical  Tubes: 

19,  33,  37,  41,  42,  45,  60,  61,  69 

Oxygen  Boiling  Outside  Vertical  Tubes: 

34,  37,  45 

Hydrogen  Boiling  Inside  Horizontal  Tubes: 

67,  74 

Nitrogen  Boiling  Inside  Horizontal  Tubes: 

43 

Hydrogen  Boiling  Outside  Horizontal  Tubes: 

48,  54 

Nitrogen  Boiling  Outside  Horizontal  Tubes: 

5,  6,  7,  18,  26,  28,  29,  36,  37,  45,  52,  54 

Oxygen  Boiling  Outside  Horizontal  Tubes: 

2,  3,  37,  45,  48,  52 

Hydrogen  Boiling  on  Horizontal  Wires: 

68,  70 

Nitrogen  Boiling  on  Horizontal  Wires: 

3,  4,  27,  28,  29,  30,  31,  32,  49,  68 

Oxygen  Boiling  on  Horizontal  Wires: 

2,  3, 4 

Hydrogen  Boiling  on  Vertical  Wires: 

73 


Nitrogen  Boiling  on  Vertical  Wires: 

3,  60,  61,  69,  71,  72 

Oxygen  Boiling  on  Vertical  Wires: 

3 

Hydrogen  Boiling  on  Inclined  Plates: 

15,  16 

Hydrogen  Boiling,  Miscellaneous  Geometry 

55 

Nitrogen  Boiling ,  Miscellaneous  Geometry 

8,  9,  10,  11,  12,  13,  14,  50,  51 

Oxygen  Boiling ,  Miscellaneous  Geometry: 

55 


OTHER  CONDITIONS 

Hydrogen  Boiling ,  Heater  Surface  Effects 

15,  16,  48,  54 

Nitrogen  Boiling ,  Heater  Surface  Effects 

18,  45,  52,  69 

Oxygen  Boiling ,  Heater  Surface  Effects: 

45,  52 

Hydrogen  Boiling,  G-field  Effect: 

1,  56,  62 

Nitrogen  Boiling,  G-field  Effect: 

8,  9,  10,  11,  12,  13,  14,  50,  51 

Nitrogen  Boiling,  Electric  Field  Effect: 

4 

Oxygen  Boiling,  Electric  Field  Effect: 

4 
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Hydrogen  Boiling,  Natural  Circulation  Effect 

57,  63 

Nitrogen  Boiling,  Natural  Circulation  Effect 

19,  38,  57,  63 

Oxygen  Boiling,  Natural  Circulation  Effect: 

38 
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APPENDIX  A 


PHYSICAL  PROPERTY  GRAPHS 


Figures  A-1  through  A-18  are  a  graphical  presentation  of  physical 
properties . 
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DENSITY  OF  OXYGEN 
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FIG.  A-4  SPECIFIC  HEAT  AT  CONSTANT  PRESSURE  OF  PARA-HYDROGEN 
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PIG.  A-5  SPECIFIC  HEAT  AT  CONSTANT  PRESSURE  OF  NITROGEN 


FIG.  A-6  SPECIFIC  HEAT  AT  CONSTANT  PRESSURE  OP  OXYGEN 


PIG.  A-7  VISCOSITY  OF  PARA-HYDROGEN 
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Thermal  Conductivity,  k,  3tu/hr-ft- 


FIG.  A-10  TrlERMAL  CONDUCTIVITY  OF  PARA-HYDROGEN 
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PIG.  A-11 
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FIG.  A-12 
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VAPOR  PRESSURE  OF  PARA-HYDROGEN 
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